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We use the cannonball (CB) model of gamma ray 
bursts (GRBs) to predict the spectral and temporal be- 
haviour of their radio afterglows (AGs). A single simple 
expression describes the AGs at all times and frequen- 
cies; its high-frequency limit reproduces the successful CB 
model predictions for optical and X-ray AGs. We analyze 
all of the observed radio AGs of GRBs with known red- 
shifts, including those of the exceptionally close-by GRB 
980425. We also study in detail the time-evolution of the 
AGs' spectral index. The agreement between theory and 
observations is excellent, even though the CB model is 
extremely frugal in the number of parameters required to 
explain the radio observations. We propose to use the scin- 
tillations in the radio AGs of GRBs to verify and measure 
the hyperluminal speed of their jetted CBs, whose appar- 
ent angular velocity is of the same order of magnitude as 
that of galactic pulsars, consistently measured directly, or 
via scintillations. 

1. Introduction 

The Cannonball Model is based on the hypothesis that 
GRBs and their afterglows are made in supernova explo- 
sions by the jetted ejection of relativistic plasmoids: "can- 
nonballs" made of ordinary baryonic matter (Dar and De 
Rujula 2000a), similar to the ones observed in quasars and 
microquasars (e.g., Mirabel and Rodriguez 1994; 1999 and 
references therein). The name cannonball (CB) originates 
in the contention that — due to a mechanism that we have 
explicitly discussed in Dado et al. 2001 — the ejected plas- 
moids stop expanding very early in the afterglow (AG) 
phase. 

The CB paradigm gives a good description of the prop- 
erties of the 7-rays in a GRB, that we modelled in sim- 
ple approximations in Dar and De Riijula 2000b. It sug- 
gests an alternative (Dar and De Rujula 2001a), which is 
rather promising (Dado et al. 2002), to the "Fe-line" in- 
terpretation of the spectral lines observed in some X-ray 



afterglows (GRB 970508: Piro et al. 1998; GRB 970828: 
Yoshida et al. 1999, 2001; GRB 991216: Piro et al. 2000; 
GRB 000214: AntoneUi et al. 2000). The model also pro- 
vides an extremely simple and successful description of 
the spectrum, and of the shape and absolute magnitude 
of the light curves of the optical and X-ray afterglows of 
all GRBs of known redshift, at all observed times (Dado 
et al. 2001, hereafter called DDD 2001). This description 
is universal, it encompasses the early optical flash of GRB 
990123, the very peculiar optical and X-ray AG of GRB 
970508, and all of the properties of GRB 980425, associ- 
ated with SN1998bw. 

In this paper we derive the CB model's predictions 
for radio afterglows, and compare them to all radio ob- 
servations in GRBs of known redshift. We also study the 
evolution of the spectral index of AGs as a function of 
time. The CB model — in parameter-thrifty and very sim- 
ple terms — passes these tests with flying colours. 

2. Summary 

In the CB model a GRB jet consists of n^^^ cannonballs, 
typically a few, each of them generating a prominent pulse 
in the 7-ray signal, as they reach the transparent outskirts 
of the shell of their associated supernova (SN). We assume 
CBs to be made of ordinary matter, mainly hydrogen, and 
to enclose a magnetic field maze, as is the case for the 
observed ejections from quasars and microquasars. The 
interstellar medium (ISM) the CBs traverse in the AG 
phase has been previously partially ionized by the GRB 
radiation and is fully ionized by Coulomb interactions as it 
enters the CB. In analogy to processes occurring in quasar 
and microquasar ejections, the ionized ISM particles are 
multiple scattered, in a "coUisionless" way, by the CBs' 
turbulent magnetic fields. In the rest system of the CB the 
ISM swept-up nuclei are isotropically re-emitted, exerting 
an inwards force on the CB's surface. This allows one to 
compute explicitly the CB's radius as a function of time 
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(DDD 2001). The radius, for typical parameters, and in 
minutes of observer's time, reaches a constant Rmax of a 
few times 10^"* cm. 

The ISM nuclei (mainly protons) that a CB scatters 
also decelerate its flight: its Lorentz factor, 7(t), is calcu- 
lable. Travelling at a large 7 and viewed at a small angle 6, 
the CB's emissions are strongly relativistically aberrant: 
in minutes of observer's time, the CBs are parsecs away 
from their source. For a constant CB radius and an ap- 
proximately constant ISM density, 7(t) has an explicit an- 
alytical expression, as discussed in Appendix I. Typically 
7 = 7(0)/2 at a distance of order 1 kpc from the source, 
and 7(0) ^ 10^. Due to a limited observational sensitiv- 
ity, GRBs have been detected only up to angles of a few 
times 1/7(0). 

The ISM electrons entering a CB are caught up and 
bounce off its enclosed magnetic domains acquiring a 
predictable power-law energy spectrum, as we argue in 
Section 3. In the CB's rest system due/dE oc below 
an energy Eb(t) ~ 7(t)moC^, steepening to duc/dE cx 
E~(P+^\ with p ~ 2.2, above this energj^. The energy 
Eb does not correspond to the conventional synchrotron 
"cooling break" but to the injection bend at the energy at 
which electrons enter the CB with a Lorentz factor 7(1). 
In Section 4 we discuss the observational support of the 
existence of the injection bend, which is strong. Given the 
very large magnetic and radiation energy densities in the 
CB, the usual cooling break (at the energy at which the 
energy-loss rate due to synchrotron emission and inverse 
Compton scattering equals that due to bremsstrahlung, 
adiabatic losses and escape) happens only at subrelativis- 
tic energies, as discussed in Appendix II. 

The magnetic energy-density in a CB (DDD 2001) is: 
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with Up the ISM baryon density (seen as 7np by the CB 
in its rest system). Thus, the magnetic field is: 
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where we have normalized fsync(i^) to a unit integral over 
all frequencies and the approximation is good to better 
than 8% precision in the range 2 < p < 2.6. Note that, 
for h' ^ ;/b, fgyn oc v^^ ^■'^^ i.e., it is independent of v^, for 
p = 2. For p ~ 2.2 , the extremely weak dependence of 
fgyn on //b in the optical and X-ray bands was neglected 
in DDD 2001. 

For the first ^ 10"^ seconds of observer's time, a CB is 
still cooling fast and emitting via thermal bremsstrahlung 
(DDD 2001), but after that the CB emissivity integrated 
over frequency is equal to the energy deposition rate of 
the ISM electrons in the CbQ 
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where no7 is the ISM electron density in the CB rest sys- 
tem and 77 is the fraction of ISM electrons that enter the 
CB and radiate there the bulk of their incident energy. In 
the early afterglow Eq. (||) must be modified to account 
for the fact that the bulk of the radio emission by the 
incoming ISM electrons is delayed by the time it takes 
them to cool down to energies much lower than their ini- 
tial one. This implies that Eq. (^) must be modified by 
a multiplicative "cumulation factor" C(i^, t), which is w 1 
at optical and X-ray wavelengths, as discussed in detail in 
Section 5. Two other factors, discussed in Sections 6 and 7, 
distinguish radio waves from higher-frequency emissions: 
attenuation by self-absorption and limb darkening; they 
introduce two extra factors A(,g[j^] and h^^{i',9^^), with 
a direction of emission relative to the CB's velocity 
vector. Normalized as in Eq. (|) and corrected by all these 
factors, the afterglow energy flux density of a CB is: 
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The pitch-angle averaged characteristic synchrotron- 
radiation frequency of electrons of energy E = Eb is 
(Rybicki and Lightman 1979): 



i^b(t) 
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where i>i^ — eB/(2 7rmoC) is the Larmor frequency in the 
CB enclosed magnetic field B. To a good approximation, 
in the CB rest system and prior to cumulation, absorption 
and limb-darkening corrections, the synchrotron radiation 
has a spectral shape: 
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^ In our previous work p was called (ip , referring to the proton 
spectrum, for which radiative losses are negligible. 



to be summed over n^-.^ for a jet with that number of 
cannonballs. This expression, for 3> fb and the second 
row set to unity, reproduces the optical and X-ray AG 
result discussed in DDD 2001. 

An observer in the GRB progenitor's rest system, view- 
ing a CB at an angle (corresponding to 0^^ in the CB's 
proper frame), sees its radiation Doppler-boosted by a fac- 
tor 5: 
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where the approximations are valid in the domain of in- 
terest for GRBs: large 7 and small 9. Since the CB is 



^ The kinetic energy of a CB is mainly lost to the ISM pro- 
tons it scatters; only a fraction < mc/mp is re-emitted by elec- 
trons, as the AG. 
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catching-upg with the radiation it emits, S is also the rel- 
ative time aberration: dtobs = dt^^g/J. The observed spec- 
tral energy density is modulated by a factor S^, two powers 
of 6 reflecting the relativistic forward coUimation of the ra- 
diation emitted in the CB's rest system. The AG spectral 
energy density Fobs seen by a cosmological observer at a 
redshift z (Dar and De Riijula, 2000a), is: 



obs 



(1 



il + z)iy S{t)t 



1 



(8) 



where Acais represents the absorption in the host galaxy 
and the Milky Way, F^^^ is as in Eq. (^), and Dl is the 
luminosity distance (we use throughout a cosmology with 
fl = 1 and flA = 0.7). In the CB model, the extinction 
in the host galaxy may be time dependent: in a day or 
so, CBs typically move to kiloparsec distances from their 
birthplace, where the extinction should have drastically 
diminished. 

In DDD 2001 we fit, in the CB model, the R-band AG 
light curves of GRBs. The fit involves five parameters per 
GRB: the overall normalization; 0: the viewing angle; 70: 
the t = value of the Lorentz factor; x^: the "decelera- 
tion" parameter of the CBs in the ISM; and the spectral 
index p. The value of p, obtained from the temporal shape 
of the afterglow, is in every case very close to the expec- 
tation p = 2.2, and — within the often large uncertainties 
induced by absorption — with the observed spectra from 
optical frequencies to X-rays (DDD 2001). 

In this paper we complete our previous work by mak- 
ing broad-band fits to the data at all available radio and 
optical frequencies. In so doing, we need to introduce just 
one new "radio" parameter: an "absorption frequency" v^, 
corresponding to unit CB opacity at a reference frequency. 
We set p = 2.2 so that the extension to a broad-band 
analysis does not involve an increase in the total number 
of parameters. We have to refer very often to the values 
of the parameters that our previous experience with the 
CB model made us choose as reference values. For conve- 
nience, these are listed in Table I. 

The predictions of the CB model, for typical parame- 
ters, are summarized in Fig. (|l|). The energy density spec- 
tra at radio to optical frequencies are shown, at various 
times after the GRB, in the upper panel. The spectral 
slopes before and soon after the peak frequency are 3/2 
and — (p — l)/2, as indicated. The spectra peak at a fre- 
quency at which self-attenuation in the CBs results in an 
opacity of 0(1). At frequencies well above the frequency 
i^b characterizing the injection bend, the spectrum steep- 
ens to a slope — p/2. In the figure's lower panel we show 
light curves at various radio frequencies. At large times 
and for v v^^ — which is the case at all frequencies in 
the example of Fig. (|^), whose parameters are close to 
those of GRB 000301c— they tend to t-2(p+i)/3^ this be- 



^ For ^7 > 1, the observer sees the back of the CB as it is 
coming forth towards her: cos^^-.^ < 1. The CB would actually 
hit or pass by the observer before its back is unveiled, were it 
not for the fact that its motion is decelerated. 



haviour being reached at earlier times, the higher the fre- 
quency. For <C J^b, the corresponding limiting behaviour 
is « t~''/3^ observable at low frequencies in the cases of 
GRBs 991216, 991208 and 000418. All of the above pre- 
dictions are robust: they do not depend on the detailed 
form of the attenuation, cumulation and limb-darkening 
factors. The early rise of the light curves does depend on 
such details, on which we shall have to invest a dispropor- 
tionate effort in Sections 5 to 7. 

The CB model provides an excellent description of the 
data, as discussed in Sections 4, 8 and 9. In the case of 
GRB 980425, for which the optical AG is dominated by 
SN1998bw, we used the parameters that fit its X-ray af- 
terglow (DDD 2001) and the GRB's fluence (Dar and De 
Riijula 2000a) to argue that they are not exceptional. The 
CB-model's description of the radio data for this GRB/SN 
pair is excellent: there is nothing peculiar about this GRB, 
nor about its associated supernova, as we discuss in detail 
in Section 9, along with the question of the angular sepa- 
ration in the sky of the SN and the associated CBs, which 
may have been, or may still be, observable. 

The apparent sky velocities of cosmological CBs are ex- 
tremely superluminal and their angular velocities happen 
to be of the same order of magnitude as those of galac- 
tic pulsars. This implies that CB velocities can possibly 
be extracted from their observed radio scintillations, as 
discussed in Section 12. 



3. The electron spectrum 

3.1. Numerical simulations 

The acceleration of charged particles by a moving CB 
is not substantially different from some of the cases al- 
ready studied in the literature, the acceleration of cosmic 
rays and electrons having attracted an enormous amount 
of attention since Fermi's first analysis in 1949 (for an 
excellent introduction, see Longair 1994). The most ef- 
ficient and thus promising mechanism is the "first-order" 
acceleration of particles by fast-moving shocks, extensively 
studied analitically and numerically since the pioneering 
works of Axford et al. (1977), Krymsky (1977), Bell (1978) 
and Blandford and Ostriker (1978). The analysis closest 
to the case at hand is that of Ballard and Heavens (1992), 
who studied acceleration by relativistic shocks, with the 
charged particles deflected by highly disordered magnetic 
fields, rather than, as it is generally assumed, by small 
irregularities in an otherwise uniform fleld. The "relativis- 
tic" and "chaotic" inputs are what make this work partic- 
ularly relevant to the case of particle acceleration by and 
within CBs. 

Ballard and Heavens study numerically, for various val- 
ues of a moving discontinuity's Lorentz factor ranging up 
to 7s = 5, the result of its collision with an isotropic en- 
semble of particles with 7p = 100. They flnd that, for 
7s = 5, the resulting particle energy distribution has a 
break (in this reference system) at 7 ~ 10 7p, at which 
point it steepens. The particles below the break have a 
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dominantly very forward motion: they are the ones which 
have been upscattered just once. Given this hint, it is easy 
to reproduce the numerical results in an analytic approxi- 
mation. In the shock's rest system, the energy of the par- 
ticles that have been scattered only once is equal to their 
incoming energy: the break in the spectrum seen in the 
simulations is a kinematical break occurring roughly at 
the injection energy. "Observed" in the system in which 
the shock is travelling at 7s = 5, this injection bend is very 
reminiscent of the familiar synchrotron-cooling "break" , 
but it has little to do with it; indeed, in the simulations of 
Ballard and Heavens (1992) cooling was entirely neglected. 

3.2. A simple analysis 

Consider the CB in its rest system and temporarily post- 
pone the discussion of cooling. The ISM electrons impinge 
on the CB in a fixed direction with a Lorentz factor equal 
to that of the CB in the GRB progenitor's rest system, 
7c = 7(t). The electrons not having "bounced back" off 
the CB's strong magnetic field, or having done it only 
once, retain the incoming energy, Eb(t) = 7(t)meC^, so 
that their energy distribution is: due/dE cx S[E — Eb(t)]. 
A very robust (i.e. detail independent) feature of the stud- 
ies of acceleration by relativistic shocks is that the parti- 
cles having bounced more than once acquire a spectrum 
due/dE (X E~P, with p = 2 in analytical approximations 
and p ~ 2.2 in numerical simulations. A few bounces are 
sufficient to attain such a spectrum. The CB is a sys- 
tem of finite transverse dimensions and the magnetic field 
contrast between its interior and its exterior is very large. 
Thus, we do not expect the same electron to bounce many 
times off the CB, as the latter catches up with it. The ac- 
celeration should occur mainly within the CB as charged 
particles bounce off its chaotically moving magnetic do- 
mains, and it should be very fast and efficient, since the 
injection is highly relativistic and there is no distinction 
between "first and second order Fermi" processes. The 
overall "source" spectrum of relativistic electrons is: 

^ ^ Ai(t)Eb(t)5[E-E„(t)] 



A2(t)e[E-Eb(t)] 



E 



Eb(t) 



with Ai and A2 of comparable magnitude and a time de- 
pendence which is that of the rate, riTrK^^^ciicj{t), at 
which electrons enter the CB. 



3.3. The spectrum of cooled electrons 

The electron energy loss by synchrotron radiation is: 



B2 



<7t C 



TT (mo c^)^ 

with /3 w 1 for the relativistic energies of interest and 
(Tt = 0.665 barn the Thomson cross-section. Let the rate 



at which fresh electrons are supplied by the ISM be called 
R. The electron source distribution of Eq. (^ "ages" by 
cooling so that: 
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At times longer than the synchrotron cooling time, the 
electron distribution tends to a time-independent duc/dE, 
obtained by equating to zero the l.h.s. of Eq. ( |ll|) and 
integrating it with the source function of Eq. (H): 
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(12) 



Admittedly, the process of acceleration that we have 
discussed is not well understood, our derivation is heuris- 
tic and Eq. ( p^ ) is not even a continuous function (the 
step function in Eq. (^) should not be so abrupt, the mag- 
netic energy in Eq. (00) should not have a fixed value). All 
we want to conclude from this exercise is that, when the 
probability of an electron to have been "kicked" only once 
is not negligible (Ai comparable to A2), the electron spec- 
trum has an injection bend at E ~ Eb(t), around which 
its spectral index changes by ~ 1 from 2 to ^ p -I- 1 . We 
choose to characterize this behaviour by the function: 



duc 
dE 



E- 



Vl + [E/Eb(t)]2(p- 



(13) 



Note how similar the injection bend is to a cooling break 
(also a spectral steepening by roughly one unit) even 
though their origins are so different. The observational ev- 
idence for an injection bend at the injection energy turns 
out to be strong, as we proceed to show. 

4. Evidence for an injection bend 

The injection bend induces the gradual transition in the 
spectral energy distribution described by Eq. (^, occur- 
ring at a "bend" frequency: 



(9) Vh 
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Hz, 



(14) 



where we have used the characteristic synchrotron fre- 
quency of Eq. ^ for the magnetic field of Eq. (||), and 
transposed the result to the observer's frame. 

For the reference CB parameters and z = 1, z/b(t = 
0) ~ 0.93 X 10^^ Hz, above the optical band. Since the 
product 7'^ 6 typically declines by more than an order of 
magnitude within a couple of days, the bend frequency in 
many GRBs crosses the optical band into the NIR during 
the early afterglow. In Fig. (j2|) we present the time de- 
pendence of ^'b(t) for 7o = 1250 and 750, characterizing 
the range of the observations, for various angles 0, z = 1, 
and the rest of the parameters at their reference values of 
Table I. The figures show that, depending on the param- 
eters, the bend frequency in the early AG may be above 
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or below the optical band, and, if it is above, it will cross 
it later. 

The bend frequency of the CB model is not the break 
frequency of the traditional fireball model. The time evo- 
lution of the former is given by Eq. and is different 
from that of the latter, which, prior to the "break" in the 
AG light-curve, can be shown to be t~^/^ (Granot and 
Sari, 2002). 

The evolution predicted by Eq.( |^) from a jy-o.5±o.i 
to a jy-i i±o i spectral behaviour is affected by extinc- 
tion. The early behaviour corresponds to times when CBs 
are not yet very far from their progenitors: extinction in 
the host galaxy may steepen the spectrum. After a day 
or more, when the CBs are further away, we do not ex- 
pect strong extinction in the host. So the prediction (after 
extinction in the Galaxy is corrected for) is an evolution 
from a behaviour close to — or steeper than — jy-o.5±o.i^ 
to a more universal i at later times. 

The predicted spectral behaviour has been observed, 
with varying degrees of significance, in the AG of several 
GRBs, listed in Table II. The first column is the bend 
frequency 1^° at t = 0, computed with Eq. (^4|) and the 
optical AG parameters of Table III (the density Up is ex- 
tracted from the measured x^o with use of Eq. ( p8[ ) and 
our reference Rmax and N^^g). For the listed GRBs the 
bend frequency is above the visible band at t = and the 
early AG measurements result in effective spectral slopes, 
/3(ti), not far from the expectation —0.5 ± 0.1, or some- 
what steeper. A few days later, the measured values, /3(t2), 
are compatible with the expectation — 1.1 ± 0.1. The sec- 
ond entry on GRB 990510 in Tabic II (Bcucrmann ct al. 
1999) requires an explanation. These authors argue that 
/3(t2) = 0.55 ±0.10, a result that assumes a strong extinc- 
tion correction in the host galaxy. But, after a day or so, 
we do not expect such an extinction. For the latest points 
measured by Beuermann et al. (1999), at day 3.85 (well 
after the bend), B-R = 0.98±0.07 and R-I = 0.49±0.06. 
Converting these results — without extinction — to a spec- 
tral slope yields /3(ti) = 1.11 ± 0.12, in agreement with 
expectation. 

The evolution from a softer to a harder spectrum 
should be a gradual change in time, rather than a sharp 
break, so that an AG's optical spectrum, if "caught" as the 
injection bend is "passing" should have an index evolving 
from —0.5 ± 0.1 to —1.1 ± 0.1 with the time dependence 
described by Eqs. (|) and In Fig. (|) we test this pre- 
diction in the case of GRB 970508, for the time-dependent 
value of the "effective" slope a ~ A [log Ft,] / A [log ly] , con- 
structed from the theoretical expectation in the same fre- 
quency intervals used by the observers. The actual pre- 
dicted fb(t) in Eq. (^ is obtained by use of the optical- 
AG fitted parameters (0, 70 and Xqo) that determine 7(t) 
and (S(t), and the density Up deduced]^ from Xqo and the ref- 
erence and Rmax- The data are gathered by Galama 



* GRB 970508 has a peculiar AG, whose CB-model interpre- 
tation requires an ISM density change at t ~ 1.2 observer's 
days (DDD 2001). 



et al. (1998a) from observations in the U, B, V, Rc and 
Ic bands (Castro-Tirado et al. 1998, Galama et al. 1998b; 
Metzger et al. 1997; Sokolov et al. 1998; Zharikov et al. 
1998), by Chary et al. (1998) for K band results, and by 
Plan et al. 1998 for the H band. 

In spite of considerable uncertainties in the spectral 
slopes deduced from observations (Galama et al. 1998a), 
the results shown in Fig. (||) are satisfactory: the observed 
crossing of the injection bend is in agreement with the the- 
oretical prediction, based on the fit in DDD 2001 to the 
overall R-band light curve from which the GRB 970508 
AG parameters have been fixed; no extra parameters have 
been fit. A couple of points in the lower panel do not agree 
with the prediction, but they do not agree with the obser- 
vations at very nearby frequencies reported in the upper 
panel, either. 

A complementary analysis to that in the previous para- 
graph is the study of an AG's optical spectrum at a fixed 
time at which the injection bend is crossing the observed 
frequency range, or is nearby. A spectral "snapshot" at 
such time should have the intermediate slope and cur- 
vature described by Eq. (^) for 1/ ^ fb. To test this 
prognosis, we compare in Fig. (^) the predicted spectral 
shape of the optical/NIR AG of GRB 000301c around 
March 4.45 UT (^ 3 days after burst) to its measured 
shape (Jensen et al. 2001). We have selected this GRB 
because its extinction correction in the galactic ISM is 
rather small: E(B - V) = 0.05 (Schlegel et al. 1998), and 
there is no evidence for significant extinction in the host 
galaxy (Jensen et al. 2001). The theoretical line in Fig. 
is given by Eq. (^ with the observer's 14, of Eq. (14) 
(i/b (1 z) = 1.75 X 10" Hz at t = 3 days, for the den- 
sity deduced from the value of Xoc of this GRB, and the 
reference values of N^^ and Rmax) ■ In the figure the the- 
ory's normalization is arbitrary but the (slightly evolving) 
slope of the theoretical curve is an absolute prediction: it 
is based on the fit in DDD 2001 to the overall R-band light 
curve and, once more, no extra parameters have been fit. 
The result is astonishingly good, even for the curvature 
which — given the figure's aspect ratio as chosen by the 
observers — is not easily visualized (a look at a slant an- 
gle helps). The late-time spectral slope deduced from the 
HST observations (Smette et al. 2001) around day 33.5 af- 
ter burst indicated a slope of ~ —1.1, again in agreement 
with our expectation. 

We conclude that the evidence is very strong for a 
spectral injection bend at the time-dependent frequency, 
Eq. (IlJ), predicted in the CB model. As illustrated in 
Fig. ffi and contrasted with data in Section 8, further 
evidence for the injection bend is provided by the fact that 
it is essential to the description of the observed broad-band 
spectra of GRB afterglows. 



5. The cumulation factor 

Three factors that are irrelevant in the optical and X-ray 
domains play a role in the description of the longer radio 
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wavelengths and the early radio AG. In this section we 
discuss the first one of them. 

Electrons that enter a CB with an injection Lorentz 
factor 7(t) are rapidly Fermi accelerated to a distribution 
that we have argued to be roughly that of Eq. (||). On a 
longer time scale, they lose energy by synchrotron radia- 
tion, and their energy distribution evolves as in Eq. (|ll|). 
Electrons with a large [7(t)] emit synchrotron ra- 

diation, with no significant time-delay, at the observer's 
optical and X-ray wavelengths. But the emission of ra- 
dio is delayed by the time it takes the electrons to "de- 
scend" to an energy at which their characteristic emission 
is in the observer's radio band. At the start of the af- 
terglow, when equilibrium conditions have not yet been 
reached, this implies a dearth of radio emission relative 
to the higher-frequency bands. This introduces a "cumu- 
lation factor" C(r/,t) in Eq. (|). 

Consider a fixed observed radio frequency Vobs- It cor- 
responds to a time changing frequency v = (l-|-z)j/obs/<5(t) 
in the CB system. The CB electrons preferentially emit- 
ting at this frequency (over an unconstrained range of 
pitch angles) are those whose Lorentz factor 7c satisfies 
the relation v ~ 0.22 7g i/l, in analogy to Eq. To esti- 
mate^ the time At it takes an electron to decelerate from 
7 ~ 7(t) to 7 = 7e, substitute the magnetic energy den- 
sity of Eq. (|^) into the electron energy loss of Eq.(p!o[) and 
integrate, to obtain 



At 



3mc 



1 



Up mp (TT c 7' 
[8.27 X 10^ s] 



"10"3cm"3- 

















(15) 



The function 7 is given by Eq. ( p?]) of Appendix I, which 
we may rewrite as: 



1 



1 

7o 



1 



t 



to = [5.14 X 10^ s] 



1 Mpc 



10^ 



n 2 



70 



(16) 



The electrons emitting the observed radio frequencies 
have 7e ~ 0(1), so that the proper CB times At and tg 
are of 0(1) year, corresponding to observer's times — 
foreshortened by a factor (1 -I- z)/S — of C(l) day. For 
optical and radio observations 70^0(7) there is no signif- 
icant delay in their emission. Moreover, the electron accu- 
mulation rate {rjir R^^x c 7 in the CB system) is orders 
of magnitude larger than the characteristic synchrotron 
cooling time E/(dE/dt) of Eq. (^, even for 7 - 10^ 
Thus, the optical and X-ray AG emission starts as soon 
as the CB is transparent to its enclosed radiation: for each 
CB, a few observer's seconds after the corresponding 7-ray 



^ We can solve Eq. (^ exactly for a given source spec- 
trum by the Mellin transform methods so familiar in Quantum 
Chromodynamics, but this would be unjustified: the accelera- 
tion mechanism is not understood well enough for the study of 
exact cooling solutions to be currently justifiable. 



pulse (DDD 2001). The radio signal, on the other hand, 
must await a time At for the cumulated electrons to cool 
down. 

The simple way to parametrize the frequency- 
dependent "cumulation effect" is to use the expression for 
the total number of electrons N(t) incorporated by the CB 
up to time t (Eq. (|40|) of Appendix I) and to positFl 



c(..,.)^^i<ij^e(.-At), 



(17) 



where the frequency dependence is via At — At(i^, t), and 
the sharp start at t > At is an artifact of our simplifi- 
cations. For optical and X-ray frequencies. At = and 
C(i^, t) = 1. In practice we find that, except for GRB 
980425 whose viewing angle is exceptionally large, one 
may also use within the radio band an approximation to 
Eq! (0): 



C(i^,t) - C(t) = 



1 - 



7(t) 
70 



1/2 



(18) 



which is also frequency independent. 
6. The attenuation factor 

At optical and X-ray frequencies the CB is transparent 
and, for the spectrum of Eq. (jj), the bulk of the radia- 
tion's energy is emitted around the bend frequency ;/b. At 
such relatively high frequencies, as illustrated in Fig. (|l]), 
absorption is unimportant. Thus, for optical and X-ray 
afterglows (DDD 2001) it suffices to know that all of the 
incoming electron's energy is reradiated, the spatial dis- 
tribution of the radiating electrons within the CB is ir- 
relevant. But in the radio, where absorption is important, 
the location of these electrons inevitably plays a role. In 
the next sections we argue that it is plausible that the ra- 
diating electrons be close to the surface "illuminated" by 
the ISM (§6.1), and that the values of the CB's plasma 
frequency (§6.2) and free-free absorption coefficient (§6.3) 
actually suggest that they may be relatively close to that 
surface. In §6.4 we deduce the final form of the attenu- 
ation factor in the CB model, characterized by a single 
parameter. 

6.1. Electron penetration 

Using numerical simulations, Achterberg al. (2001) have 
shown that for simple geometries the bulk of highly rela- 
tivistic particles encountering a coUisionless shock escape 
before they undergo diffusive shock acceleration. In real- 
ity, the geometry of the CB, its density distribution and 
its magnetic field distribution are very complicated, mak- 
ing the fraction of the ISM electrons that penetrate inside 
the CB, and their distribution there, very uncertain. 



The rapid onset of the radio signals from SNe is not under- 
stood, see e.g. Weiler et al. 2001. Perhaps electron cumulation 
also plays a role there. 
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Several length scales play a role in discussing the fate 
of an electron that enters the CB with 7c = 7(t)- The 
Larmor radius is Rl = moC^7/(eB), which is indepen- 
dent of 7 for B scaling as in Eq. (||). For our reference pa- 
rameters, Rl ~ 6 km is many orders of magnitude smaller 
than the CB's radius and does not play a crucial role. The 
length of an electron's curled-up trajectory as it radia- 
tively loses energy is cE/(dE/dt) or ^ 2.6 x 10^^ cm for 
the cooling rate of Eq. ( p^ and an initial 7 = 10"^. This 
is only an order of magnitude larger than the reference 
CB's radius Rmax- We have no way to estimate the typi- 
cal coherence size of a CB's magnetic domain Lb, but the 
depth D ~ {cTj LbY^^ to which an electron penetrates, 
even for a relatively simple magnetic mess (Lb not much 
smaller than Rmax) is smaller than the CB's radius. For 
Lb as small as Rl, D ~ 4 x 10^" cm, some four orders of 
magnitude smaller than Rmax- Even this concrete value is 
uncertain, for it depends on the surface magnetic field as 
B~'^/*, and the surface B- value may be different from that 
of Eq. which is a volume average. 

In addition to all of the above uncertainties, it is possi- 
ble that a CB's illuminated working surface be turbulent, 
and harbour fast plasma motions, if only to establish local 
charge neutrality, which is disrupted as electrons and pro- 
tons penetrate the CB to different depths. We conclude 
that the the fraction of ISM electrons that enter inside 
the CB may be small and the synchrotron-radiating elec- 
trons may be concentrated close to the CB's surface, as 
opposed to be acquiring a uniform distribution over the 
CB's volume. 



6.2. The plasma frequency 

The plasma frequency in a CB with an average free inner 

free 



electron density n^'''^'' is 



1/2 



28 



^frec 



10^ cm-3 



1/2 



MHz. 



(19) 



For a fully ionized CB 11^'°° = fie , to whose reference value 
we normalized the above result (the fraction of electrons 
swept up from the ISM is small, relative to the total num- 
ber in the CB -to whose free fraction Eq. ( [T9| ) refers- but 
the CB is highly ionized, as shown in Appendix III). 

For v < Vp the radio emission is completely damped 
within a typical length ~ c/[i^p — v^]^^^, much smaller 
than the CB's radius. At very early times, 5 ~ 10"^ and 
the Doppler-boosted value of i^p falls in the low end of the 
observed range of radio signals, where a sharp cutoff is not 
observed. We must conclude that the (small fraction of) 
radiating electrons is located in a CB surface layer whose 
total electron density (dominated by the thermal electron 
constituency) is smaller than our reference average value, 
a one order of magnitude reduction being comfortably suf- 
ficient to move the value of Vp to a position below the cur- 
rently observed frequencies. We have explicitly checked 
that our fits do not improve significantly with the inclu- 
sion of Vp as a free parameter: the minimization procedure 



always "gets rid" of the fit by choosing it somewhat be- 
low the reference value of Eq. ( [l^ ) , and below the lowest 
measured frequencies. 

6.3. Free- free attenuation 

At the MHz frequencies in the CB system corresponding 
to the observed radio frequencies, the synchrotron emis- 
sion is strongly attenuated by free- free absorption (inverse 
bremsstrahlung) in the CB; free-free absorption dominates 
over self-synchrotron absorption, as shown in Appendix 4. 
For a hydrogenic plasma, the free-free absorption coeffi- 
cient at radio frequencies is: 



0.018 gff ho hi T"3/2 j,-2 ^^^-1 



(20) 



where lii ~ ho is the free ion density in the CB, in units of 
cm^'^, T is the plasma temperature in degrees Kelvin, v 
is in Hertz and gs is a Gaunt factor for free-free emission, 
of 0(10) at the relevant frequencies. 

The opacity t,^ of a surface layer of depth D is: 



.(D,t) 



dr = x^[t]D. 



(21) 



Equilibrium between photoionization of atomic hydrogen 
in the CB by synchrotron radiation and recombination 
of free electrons and protons to hydrogen keeps the CB 
plasma partially ionized during the observed AG. The 
Coulomb relaxation rate in the CB is very fast because 
of its high density. Consequently, the CB plasma is ap- 
proximately in quasi thermal equilibrium. Because of the 
exponential dependence of the Saha equation on temper- 
ature, and the high ionization rate, the CB's temperature 
is kept practically constant around a few eV, and the ion 
density and free electron density become proportional to 
7(t) (Appendix III). Using Tq = 10^ K and the reference 

lO'^cm-3 in Eqs. (|0) and 



average densities Uo 
(^, we obtain: 

T^(D,t) ~ 1.4 X 10^ 



Ui 



D 



Rn 



1 GHz 



-1 -3/2 



105 



(22) 



which is very large for D ^ ©(Rmax)- A reduction in 
surface or average CB density of one order of magnitude 
or more — which, as we have seen, renders unobservable 
the unobserved plasma-frequency cutoff — reduces t^, by 
two orders of magnitude, or more. For D Rmax/ 10, 
this would make t(D), as required, of order unity at the 
peak frequency ~ 10^ GHz of the early-time spectrum of 
Fig. (j^) (at which time the observed and CB frequencies 
differ by a factor 6/{l + z) ^ 10^). 

The conclusion is that a reasonable deviation of the 
properties of the CB from their reference bulk average 
values (a reduction of the total number-density of free 
electrons in a synchrotron-emitting surface layer) implies, 
not only that the plasma-frequency break is not observable 
in the current data, but also that the magnitude of the 
free-free attenuation is the required one. Our ignorance of 
the depth, temperature and density of ions and electrons 
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in the radio-emitting surface of a CB can be absorbed into 
a single parameter: a characteristic absorption frequency, 
i/a, in the opacity of Eqs. (|o[ |l[ |2|): 

(23) 



The frequency dependence of the free-free attenuation, 
Xi/ c)c is fairly well supported by the observed ra- 

dio spectra at their lowest frequencies, as our comparisons 
with observations in Sections 8 and 9 demonstrate. 



6.4. Attenuation in slabs and spheres 

We do not know a priori the geometry of the working sur- 
face from which a CB's synchrotron radiation is emitted. 
In the case of optical AGs this is immaterial, for the CB is 
transparent to radiation at the corresponding CB-system 
wavelengths: the bulk of the radiation energy is emitted at 
these frequencies. For the case of radio AGs, attenuation 
is important and the shape of the emitting surface layers 
plays a role: the expression for attenuation as a function 
of opacity is geometry-dependent. 

For a planar-slab geometry, the familiar expression for 
the attenuation is: 



AM = 



1 - c-^- 



(24) 



For the emission from a sphere of constant properties, we 
obtain: 



AH = 



3t2-6t^ + 6[1-^ 



(25) 



while for the emission from a thin spherical surface, the 
result of Eq. (|2^ ) is recovered. 

For the sake of definiteness, we adhere to CBs that are 
spherical in their rest system. This means that, as the fre- 
quencies increase and the CB evolves from being opaque 
to being transparent, we should use an attenuation evolv- 
ing from Eq. ( p5| ) to Eq. (p^. Rather than doing that, we 
have checked explicitly that our results are insensitive to 
the use of one or the other form, and used the simpler one. 

7. The illumination and limb-darkening factor 

Consider a spherical CB in its rest system. It is "illumi- 
nated" by incoming ISM electrons only in its "front" hemi- 
sphere. If observed at an angle 6^,^ 7^ 0, a fraction of the 
"dark" CB is also exposed to the observer, like the Moon 
in phases other than totality. For radio waves — to which 
the CB is not transparent — these geometrical facts play 
a non-trivial role. 

Place the direction of the CB motion, or of its illu- 
mination, at {OjCp) — (0,0); at a direction = (0,0,1) 
in Cartesian coordinates. The normal to a sphere's sur- 
face point at {9,(p) is rig ~ (cos 6* sin sin6', cos cos 0). 
The observer is in the direction (0,6^^), where we have 
taken the liberty to label what in this parametriza- 
tion is an azimuthal angle; the corresponding unit vector is 



— (sin 0(^3 , 0, cos 0(,g). The relation between 9^^ and 
the terrestrial observer's viewing angle is that of Eq. (|^). 

When attenuation plays a significant role, an element 
of a CB's surface reemits an amount of energy propor- 
tional to the cosine of the illumination angle: rii ■ ris. 
Because of the limb-darkening effect, the reemitted ra- 
diation depends on the cosine of the angle between the 
surface element and the observer: rig ■ n^B = A simple 
characterization of the functional form of the limb dark- 
ening effect (see e.g. Shu, 1991) is: 



(26) 



The combined effect of illumination and limb darkening is 
an emitted radiation proportional to: 

E(cos6'^b) — I dcosO i d0F(ns • n^g) rii • rig 

= i[2-f (2 + ^-0^3) cos0^3+sin0^3]. (27) 



An excellent and simple approximation to Eq. ( |27| ) is: 
E(x) = E(l)^[4-t-x][l-t-x], (28) 

with X = cos 0(^3 . 

For negligible self- attenuation A[i^] = 1, as in the op- 
tical, there is no limb darkening and illumination effect. 
As absorption becomes increasingly important for longer 
wavelengths, the effect becomes fully relevant. We inter- 
polate between these two extremes by writing: 

L,3(^.,cos0,3) ^ AH + (1 - AH) (29) 

to obtain the overall illumination and limb-darkening cor- 
rection factor to the energy flux density of Eq. (0). 



8. Broad band spectra: radio and optical results 

In practice, it is not an effortless task to test a prediction 
for an AG's spectrum extending, as in the upper panel of 
Fig. (|l|) , to all measured wavelengths from radio to X-rays. 
The problem is not related to the model, but to the data. 
First, the corrections due to absorption, particularly in 
the host galaxy, are frequency-dependent and notoriously 
difficult to ascertain with confidence. Second, the integra- 
tion times employed in the radio observations are long, so 
that the theoretical prediction varies within the time win- 
dow, and so do the optical energy flux densities, measured 
over much shorter periods, as well as some of the radio 
observations themselves. Unavoidably, this will make our 
spectral flgures look a bit peculiar, with two theoretical 
curves bracketing the expectations, and various observa- 
tional points at the same frequency. 

We study the AG light-curves and broad-band spectra 
of all GRBs with known-redshift whose AG was measured 
both in the radio and optical bandsj^. 

^ The domain extending from the optical to the X-ray regime 
— but for the early injection bend discussed in Section 4 — is 
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Our predictions are given by Eq. (||), fit to the opti- 
cal and radio observations. The fitted parameters are the 
overall normalization, 7(0), 0, the deceleration parameter 

(whose meaning and role are reviewed in Appendix I) 
and the CB self absorption frequency of Eq. (|2^) . We 
found in DDD 2001 that p is very narrowly distributed 
around its theoretical value p = 2.2, and we fix it to that 
value for all GRBs in the current analysis. Thus, the total 
number of parameters in our broad-band fits is the same 
as we used in DDD 2001 to describe just the R-band light- 
curve. 

The values of the parameters, listed in Table III, are 
very similar to those deduced in DDD 2001 by fitting only 
the R-band optical data with the high-t/ limit of Eq. (^. 
The small differences are due not only to the use of ra- 
dio data and optical bands other than R, but also to the 
inclusion of the effects of the injection bend in the CB 
synchrotron AG global formula, Eq. (||), and (to a small 
extent) to the use of a fixed p = 2.2. The results show 
that the theory agrees with observations both at radio 
and optical wavelengths. For some GRBs a slightly bet- 
ter fit to the radio data is obtained if the absorption fre- 
quency Va. is best fitted to the radio data alone or if a fitted 
power-law dependence on time is used for the CB opac- 
ity instead of Eq. (p3|), with all other parameters taken 
from the global fit. Because of scintillations, and of the 
very detail-dependent character of our prediction for the 
time dependence of a CB's opacity, it is difficult to as- 
sess whether or not the slightly improved values are 
significant or not. 

Notice in Table III that the distributions of parameters 
are fairly narrow, in particular for 70 . Of particular inter- 
est, since it can be predicted, is the distribution in 7 x 0. 
Since AGs are discovered at optical and X-ray frequencies, 
the angular distribution is that of the high frequency limit 
of Eq. (|). For smaU 6*, dN/d6l cx / . This dis- 
tribution has a maximum at 76 ~ 0.37 and a median at 
^ 0.5. In Table III there are four cases with 76* below the 
median and five above. The worst "outlier" in 6 is much 
less so in 76*. The conclusion that this distribution is per- 
fectly compatible with the expectation can also be reached 
from Fig. (39) of DDD 2001, whose results were obtained 
from only optical data, but for which the statistics is a bit 
better. 

We first discuss the broad-band spectra and light 
curves of three representative GRBs: 000301c, 000926 and 
991216. The optical AG of GRB 000301c is practically 
unextinct, that of GRB 000926 has strong extinction in 
the host galaxy (e.g., Fynbo et al. 2001) and that of 
GRB 991216 has strong extinction both in the host galaxy 
and in ours (e.g., Halpern et al. 2000). We discuss GRBs 
991208, 000418, 000510, 990123 and 970508 in slightly less 
detail. The apparently special case of GRB 980425 is dis- 
cussed separately in the next chapter. 



compatible with the expected behaviour v ~ v ^'^ (DDD 
2001). 



8.1. GRB 000301c 

For this GRB we fit the radio data of Berger et al. (2000) 
and the optical data of Garnavich et al. (2000b), Jensen 
et al. (2001), Masetti et al., (2000), Rhoads and Fruchter 
(2001) and Sagar et al. (2001). Our results for the light 
curves at all observed optical and radio frequencies are 
gathered in Fig. (||), which is representative of the trends 
seen in all GRBs. The narrowly spaced lines in the figure 
are the optical light curves for — from top to bottom — 
the K, J, I, R, V, B and U bands. Their very satisfactory 
comparison with data is reported in Fig. (|^). The results 
for the radio AG are the more spaced lines in Fig. (||), 
which correspond — from top to bottom at the figure's 
left side— to frequencies of 1.43, 4.86, 8.46, 15, 22.5, 100, 
250 and 350 GHz. Their very satisfactory comparison with 
observations is reported in Figs. (|^) to (p^. Notice that 
all features of the data have precisely the trends summa- 
rized in Fig. (H). In Figs, ( p]) and ( |l2| ) we present the 
complementary information, by comparing our fits with 
the observations for the radio-to-optical spectra of GRB 
000301c in four radio time- integration brackets; 1 to 5, 
5 to 10, 10 to 20, and 20 to 30 days. The pronounced 
peaks are at (observer's) frequencies for which the opacity 
of Eq. (|23| ) is ~ C'(l)- The injection bend at a higher 
frequency is clearly visible, it is responsible for the agree- 
ment between the radio and optical magnitudes and fre- 
quency trends. The two curves in these figures, and many 
later ones, refer to the expectation at the two times which 
bracket the actual radio observation. The results are quite 
satisfactory. 

8.2. GRB 000926 

We have made a global fit to the NIR/optical data (Di 
Paola et al. 2000; Fynbo et al. 2001, Harrison et al. 2001; 
Price et al. 2001; Sagar et al. 2001) and the radio data 
(Harrison et al. 2001) on this GRB. In Fig. (|l|) we com- 
pare the fitted CB-model predictions with the measured 
light curves for the I, R, V, B and U bands, after subtrac- 
tion of the host galaxy and SN contributions (DDD 2001). 
The theoretical predictions were corrected for galactic ex- 
tinction E(B - V)=0.0235 (Schlegel et al. 1998) and for the 
estimated extinction in the host galaxy, E(B - V)=0.40 
(Harrison et al. 2001). In Figs. (|l^ ) to (|l^ ) we present the 
radio light curves for six frequencies ranging from 98.48 
to 1.43 GHz. 

In Figs. (|l^) to (|l^) we make the complementary com- 
parison of theory and observations for the radio-to-optical 
spectra, in six time intervals extending from 0.8 to 100 
days. The results, in spite of the crude estimate of extinc- 
tion in the host galaxy and the scintillations so clearly 
visible in the radio light curves, are satisfactory. 

8.3. GRB 991216 

The NIR/optical data for this GRB are from Halpern et 
al. (2000) and Garnavich et al. (2000a); the radio data 
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from Frail et al. (2000b). In Fig ( |20| ) we present the com- 
parison between the measured hght curves for the K, J, 
I, R bands, after subtraction of the host galaxy and SN 
contributions (DDD 2001), and the fitted CB model pre- 
dictions. The predictions were corrected for extinction in 
the host galaxy and ours, as estimated by Halpern et al. 
(2000): E(B - V)=0.40. In Figs. (|l|) to (|2|) we present 
the radio light curves at six frequencies from 350 to 1.43 
GHz. In Figs. (|2^) to (|2^) we make the complementary 
comparison of theory and observations for the radio to 
optical spectra, in six time intervals extending from 0.44 
to 80 days. The results are once again satisfactory. 

8.4. GRB 991208 

We fit the NIR/optical data (Castro-Tirado et al. 2001; 
Sagar et al. 2000) and the radio data (Galama et al. 2000) 
on the AG of GRB 991208. In Fig (H) we present the com- 
parison between the measured light curves for the I, R, V 
and B bands, and the fitted CB model predictions, after 
subtraction of the host galaxy and SN contributions (DDD 
2001). The theoretical predictions were corrected only for 
the small galactic extinction E(B - V)=0.016 (Schlegel et 
al. 1998) in the direction of this GRB, there being no spec- 
tral evidence for optical extinction in the host galaxy. In 
Figs. ( p7| ) to (30) we also present the radio light curves at 
100, 86.14, 30, 22.49, 14.97, 8.46, 4.86 and 1.43 GHz. In 
Figs. ( ^ ) and ( ^2| ) we make the complementary compari- 
son for the radio to optical spectra in three time intervals 
extending from 2 to 14.3 days. The results are satisfactory. 

8.5. GRB 000418 

The NIR/optical data are from Klose et al. (2000) and the 
radio data from Berger et al. (2001a). In Fig. ( |33| ) we com- 
pare the fitted CB-model predictions with the measured 
light curves for the R-band, after subtraction of the host 
galaxy and SN contribution (DDD 2001). The theoretical 
predictions were corrected for galactic extinction and for 
extinction in the host galaxy as estimated by Berger et 
al. (2001a): E(B - V)=0.40. In Figs. (||) and (H) we also 
present the radio hght curves at 22.5, 15, 8.46 and 4.86 
GHz. In Fig. ( |3^ ) we make the complementary compari- 
son for the radio to optical spectra, in two time intervals 
extending from 9.5 to 100 days. The results are satisfac- 
tory. 

8.6. GRB 990510 

The NIR/optical data were gathered by Beuermann et al. 
(1999), Harrison et al. (1999) and Stanek et al. (1999) and 
the radio data by Harrison et al. (1999). In Fig. (^^ we 
present the comparison between the measured light curves 
for the I, R, V, B bands, after subtraction of the host 
galaxy and SN contribution (DDD 2001), and the fitted 
CB model predictions, corrected for Galactic extinction 
(E(B - V) =0.203, Schlegel et al. 1998) and for extinction 



in the host galaxy as estimated by Stanek et al. (1999). 
In Figs. (^) and (39) we present the radio light curves 
at 13.7, 8.6 and 4.8 GHz. In Figs. (||) and (^ we also 
make the complementary comparison of theory and obser- 
vations for the radio to optical spectra three time intervals 
extending from 1 to 40 days. The agreement between the- 
ory and observations is very good although its significance 
is limited by the sparse radio data. 

8. 7. GRB 990123 

We have fit the NIR/optical data (Castro Tirado 1999; 
Fruchter et al. 1999; Galama et al. 1999, Holland et al. 
2000; Kulkarni et al. 1999a) and the radio data (Galama et 
al. 1999; Kulkarni et al. 1999b) for this GRB. In Fig. ^ 
we present the comparison between the fitted CB model 
predictions — assuming a constant ISM density after 0.1 
observer's days^ and after subtraction of the host galaxy 
and SN contributions — with the measured light curves 
for the K, I, R, V, B and U bands. The theoretical predic- 
tions were corrected for the small Galactic extinction in 
the GRB direction (E(B - V)=0.016, Schlegel et al. 1998) 
but not for extinction in the host galaxy, since there is 
no spectral evidence for significant extinction there. In 
Fig. (^) we present the radio light curves at 15 and 8.46 
GHz. In Figs. ( |43| ) to ( ^ ) we make the complementary 
comparison of theory and observations for the radio to op- 
tical spectra, in four time intervals extending from 0.1 to 
20 days. The agreement between theory and observations 
is good despite the limited available data on the radio AG 
and its modulation by scintillations. 

8.8. GRB 970508 

The optical (and X-ray) AG of GRB 970508 is the only 
one so far that has been seen to rise and fall very sig- 
nificantly (e.g., Garcia et al. 1998; Galama et al. 1998b; 
Pedersen et al. 1998; Schaefer et al. 1997; Sokolov et al. 
1998; Zharikov et al. 1998). In DDD 2001 we have shown 
that a CB model fit to this AG fails, if one assumes — like 
in all our other fits — a constant ISM density. However, we 
have argued there that GRB progenitors are presumably 
located in super-bubbles of 0.1 to 0.5 kpc size. There may 
be instances in which the jet of CBs, after travelling for 
such a distance, does not continue onwards to a similarly 
low-density halo region, but encounters a higher-density 
domain. Indeed, we have shown that a fairly satisfactory 
fit to the optical (and X-ray) AG is obtained upon assum- 
ing an upwards jump in density by a factor 2.2 at t ~ 1.1 
day after burst. This jump occurs before the first available 
data points on the radio AG (Galama et al. 1998a; Frail 
et al. 2000b). Therefore, we have fitted the optical data 
and the radio data with the ISM density profile that was 
fitted to the R-band light curve. 

In Fig. ( ^5| ) we present the comparison between the 
measured light curve for the I,R,V and B bands after sub- 



the earlier optical data are discussed in DDD 2001. 
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traction of the host galaxy and SN contribution (DDD 
2001). The theoretical predictions were corrected for the 
small galactic extinction in the GRB direction (E(B - 
V)=0.016, Schlegel et al. 1998) but not for extinction in 
the host galaxy, since there is no spectral evidence for sig- 
nificant extinction there. In Figs. (46) and ( ^7| ) we also 
present the radio light curves at 8.46, 4.86 and 1.43 GHz. 
In Figs. ( ^7| ) to (^) we make the complementary com- 
parison of theory and observations for the radio to optical 
spectra, in five time intervals extending from 0.12 to 470 
days. The results are quite satisfactory. 

8.9. Commentary 

In DDD 2001 we demonstrated that, in the CB model, 
the spectral index in the optical to X-ray domain could 
be extracted from the time-dependence of the optical light 
curves. The fits resulted in a = p/2 ~ 1.1 for all GRBs of 
known redshift. This result is in good agreement with the 
observed late spectral observations. We have learned in 
this section that the CB model also provides an excellent 
description of the AG spectra in the broader band that 
includes the radio data. Only one new parameter, z^a, is 
involved in the extension to the broader band. And this 
fitted parameter and the injection bend — at its predicted 
frequency and time-dependent position — bring about the 
agreement between the different magnitudes and spectral 
trends of the radio and optical domains. 

In some of our fits to broad band spectra, such as 
the earhest data on GRBs 000301c, 991216 and 990123 
in the upper panels of Figs. ([TTI), ( p3| ) and (|3|), respec- 
tively, the theoretical curve is an underestimate of the 
low-energy spectral intensity. In other cases, such as GRBs 
000926, 991208 and 980425, the spectral fits are excellent 
at all times. The lowest frequencies and earliest times are 
the most dependent on our simplifications concerning the 
GRB geometry, density profile, self-absorption, cumula- 
tion and limb- darkening. We would have been surprised 
if these simplifications worked even better than they do, 
and the fits do improve if we remove our approximation of 
a fixed spectral index p = 2.2. But our aim in this paper 
is not to obtain spectacularly good fits, but to demon- 
strate that, even in the simplest approximations, the CB 
model provides a good description of the broad-band data. 
The analysis of the lowest radio frequencies at the earliest 
times brings forth a plethora of details that are not of fun- 
damental interest: our ultimate goal is not to understand 
these details, but to investigate what the origin of GRBs 
actually is. 

9. SN1998bw and GRB 980425 

The time and position of the peculiar gamma ray burst 
980425 (Soffita et al. 1998) coincided with supernova 
SN1999bw (Tinney et al. 1998) in the spiral galaxy ESQ 
184-G82, at a nearby z = 0.0085 (Tinney et al. 1998; 
Sadler et al. 1998; Galama et al. 1998c; Lidman et al. 1998; 
Iwamoto et al. 1998). Iwamoto et al. (1998) estimated the 



the core collapse of SN1998bw to have happened within —2 
to +7 days of GRB 980425. The BeppoSAX Narrow Field 
Instrument (NFI) located lOh after burst an X-ray source 
coincident in position with SN1999bw that declined slowly 
with time between April and November 1998 (Plan et al. 
2000). A posteriori statistics indicate a very low chance 
probability (< 10~^) of a GRB being so nearly coincident 
in position. But despite how close — if it was associated 
with SN1999bw— the progenitor of GRB980425 was to us 
(38 Mpc for H = 65kmMpc~^ the 7-ray fiuence indi- 
cated only 7 X 10^* erg equivalent spherical energy release 
in 7 rays, much smaller than ~ 3x 10^"^ erg, the mean value 
for the score of other GRBs with known — cosmological — 
redshifts. 



9.1. SN1998bw: the accepted lore 

Like its accompanying GRB, SN1998bw was also claimed 
to be a very peculiar radio supernova (e.g. Kulkarni et 
al. 1998). Over the past twenty years approximately two 
dozen SNe have been detected in the radio: 2 Type lb, 5 
Type Ic, and the rest Type II. A much larger list of more 
than 100 additional SNe have low radio upper limits (for a 
review see, e.g., Weiler et al. 2000 and references therein). 
Type Ib/c SNe are fairly homogeneous in their radio prop- 
erties, but SN1998bw0 had a peak 6-cm radio luminosity 
of 8 X 10^^ ergs Hz "\ that is 20 to 40 times brighter 
than other radio Type Ib/c SNe, which fall typically in 
the range 1.4 - 2.6 x 10^* ergs"^ Hz'^ SN1998bw also 
reached a high radio luminosity earlier than any known 
SN. Simple arguments based on the brightness tempera- 
ture of its radio luminosity (e.g., Readhead 1994) required 
the radiosphere of SN1998bw to have expanded surpris- 
ingly fast, at > 200, OOOkms"^, at least during the first 
few days. Its unusually high optical and radio luminosities 
and its extraordinarily large initial speed of expansion led 
many authors to conclude that SN1999bw was a hyper- 
nova (Paczynski 1998) rather than a peculiar supernova 
(e.g., Iwamoto et al. 1998). 

9.2. The pair SN1998bw/GRB 980425 in the CB 
model 

In Dar and De Riijula (2000a) we argued that the only 
peculiarity of SN1998bw was that it was viewed very near 
its axis. The peculiarity of GRB 980425 was its nearness, 
that allowed for its detection at an angle, ~ 8/7(0) 
unusually large relative to the other GRBs of known red- 
shift, for which 6 ^ 1/70- These facts conspired to produce 
a "normal" GRB flucnce, and resulted in an optical AG 
dominated by the SN. In DDD 2001, we demonstrated 
that the X-ray AG of this GRB was also "normal" : it has 
precisely the light curve (in shape and normalization) ex- 



^ SN1998bw was classified initially as Type lb (Sadler et al. 
1998), then Type Ic (Patat and Piemonte 1998), then peculiar 
Type Ic (Kay et al. 1998), then, at an age of 300 to 400 days, 
again as Type lb (Patat et al. 1998). 
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pected in the CB model if the X rays are produced by the 
CBs and not, as the observers assume (Pian et al. 2000), 
by the supernova. 

In the CB model (Dar and De Riijula 2000a), the 
gamma-ray fluence of GRBs at large viewing angle 



hio' » 1) 



IS oc 



cx 



The radio AG spectral en- 



ergy density is proportional to as implied by 

Eqs. (^, H, H), the dependence Vh « 7^ and the relation 
i^obs oc 6. As a function of time, the AG peaks when 
76* ~ 1, so that 7 ~ (5 ~ 1/0 and the peak value is pro- 
portional to B^"^. Because its proximity and large viewing 
angle "conspired" to make GRB 980425 appear "normal" 
in gamma rays, its peak radio intensity should have been 
enhanced by a factor ~ (6'/mrad)^ relative to that of ordi- 
nary GRBs. Thus, for 9 ^ 8.3 mrad, as estimated for GRB 
980425 in Dar and De Riijula 2001, its expected peak ra- 
dio intensity is ^ 60 times larger than that of ordinary 
GRBs. Observationally, it is 50 to 100 times larger. 

In Figs. ( |50| ) to (??) we show our CB model fits to the 
temporal and spectral behaviour of the radio afterglow of 
GRB 980425. The fit parameters (in particular the large 
observation angle 6) are quite close to the ones that ex- 
plain its GRB fiuence (Dar and De Riijula 2000a), and 
its X-ray afterglow (DDD 2001). These figures show how, 
in the CB model, the radio AG of GRB 980425 also has 
a "normal" magnitude and shape. That is, once more, if 
the radio AG is produced by the CBs and not by the SN, 
unlike, once again, it is generally assumed (e.g., Kulkarni 
et al. 1998; Li and Chevalier; Weiler et al. 2000). 

In the case of GRB 980425 the relatively large viewing 
angle and the subsequently small Doppler factor imply 
that, at late times, even the radio frequencies are above 
the injection bend. The large v behaviour in Figs. (p3[) 



and {l 



IS 



J.-P/2 



Also, the late time trend of the 



radio light curves in Figs. (|50| ) and (51) approaches the 
asymptotic t^fp+i)/^ ~ t^^.i^ 

For GRB 980425 the radio data are overwhelmingly 
more abundant and precise than the X-ray data, and it is 
interesting to check what the prediction for the X-ray light 
curve is, if the input parameters are those determined in 
the radio fits. This is done in Fig. (^5|) for two values of 
the electron spectral index p. For our fixed choice, p = 2.2, 
the prediction misses the data by a factor ~ 20. There are 
two excuses for that. First, since GRB 980425 is seen much 
more "sideways" than other GRBs, and its Doppler factor 
S is much smaller than usual, the cumulation, illumination 
and limb-darkening factors play a bigger role than usual. 
These factors involve many simplifying assumptions (such 
as spherical symmetry) and significantly affect the nor- 
malization of the radio AG, but not that of the X-rays. 
Second, the extrapolation from radio to X-rays is over 
some 10 orders of magnitude in frequency, and a small 
change in the spectral photon's slope, (p — l)/2, entails a 
very large change in relative magnitude, as can be seen in 
Fig. ( |55| ) by comparing the p = 2.2 and p = 2 curves. 

The X-ray light-curve of GRB 980245 is essentially flat 
in the time-interval of the first four observed points (Pian 
et al. 2000), while the corresponding data for all other 



GRBs fall with time much faster. The last observational 
point in Fig. ( ^sf ) , a preliminary result from XXM Newton 
(Pian 2002) and Chandra (Kouveliotou 2002), falls pre- 
cisely in the expected subsequent fast decline (predicted 
in DDD 2001) and definitely not in a naive power-law 
extrapolation. The peculiar light curve is a consequence 
of the large observing angle (Dar and De Riijula 2000a). 
For the reasons stated in this paragraph and the preced- 
ing one, we consider the prediction of the X-ray fluence 
completely satisfactory. 

In DDD 2001, on the basis of the very meager X-ray 
data, we argued that the last optically-measured point of 
the SN1998bw/GRB 980425 pair, at day 778 (Fynbo et 
al., 2000), was due to the CB's AG and not to the super- 
nova. Redoing the analysis with the input of the abun- 
dant radio data, we must now revise this conclusion. In 
Fig. ( p4[ ) we show the result, with inclusion of the late op- 
tical measurement. This point lies more than two orders of 
magnitude above the predicted CB's AG: it must be due 
to the SN. We do not have an explanation — specific to the 
CB-model — of the fact that this point also lies somewhat 
above the expectation based on '"^^Co decay (SoUerman et 
al. 2000). 

We are claiming that long duration GRBs are asso- 
ciated with a good fraction of core-collapse SNe. Yet, 
SN1998bw was one of the brightest in its class. The appar- 
ent contradiction may be dispelled by the increasing evi- 
dence that SN explosions are fairly asymmetric. It is quite 
conceivable that, viewed very close to their "CB axis" SNe 
appear to be brighter than when observed from other di- 
rections. 

The conclusion is twofold. GRB 980425 is, in every 
respect, normal (z and 6 being chance variables). And, 
deprived of very abnormal X-ray and radio outputs — 
which are not due to the supernova, but to its ancillary 
GRB— SN1998bw loses most of its "peculiarity" . 

9.3. Superluminal motion in SN1998bw/GRB 980425 

The transverse projected velocity in the sky of a CB rel- 
ative to its parent SN is, for large 7 and small 9: 



(30) 



which, for typical parameters, is extremely superluminal. 
The resulting angular separation at time t is: 



Aa(t) 



1 



VcB(t')dt' 



(31) 



where Da = (1 + z)^ Dl is the angular distance to 
the SN/CB system. In Fig. (|5|) we show A(t) for 
SN1998bw/GRB 980425 with our parameters fit to the 
corresponding radio data (for our adopted Hq, Dl — Da — 
38 Mpc). In Dar and De Riijula (2000a) we argued that 
this separation was sufficient to justify a dedicated effort 
to search for a "binary" source. It is interesting to discuss 
what the situation is with the data currently available. 
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The most accurate determination of the position in 
the sky of the SN1998bw/GRB 980425 system is based on 
the radio observations made with the Austrahan Telescope 
Compact Array (ATCA, Wieringa et al. 1998). RecaU that 
in the CB model the radio coordinates are those of the 
CB (GRB 980425 was a single-pulse GRB, that is, it had 
a single dominant CB). In days 3, 4 and 10 the source 
is reported to be at (RA 19:35:03.31, Dec -52:50:44.7). 
In the subsequent 33 observations, ranging from day 12 to 
day 790, the position is (RA 19:35:03.32, Dec -52:50:44.8), 
some 0. 18 away from the original determination, but not 
inconsistent with the observational uncertainty of 0. 1. In 
the penultimate observation at day 320 the source has 
faded to the point that it is not observable in 2 out of 6 
frequencies, and in the last observation at date 790 there is 
no clear sighting at any frequency. The predicted values of 
Aa from Eq. (30) at some relevant dates are 12, 158, 183 
and 292 mas at days 12, 249, 320 and 790, respectively. 
These results, the observational error, and the fact that 
the ATCA observers were not trying to follow the source's 
motion imply that their results are insufficient to claim 
either that the early change of position was significant, or 
that a motion of the CB comparable to the predicted one 
is excluded. 

Observations of the vicinity of the source of GRB 
980425 were made with the Hubble Space Telescope (HST) 
at day 778 (Fynbo et al., 2000), with a tiny astrometric 
uncertainty of 0. 018, and pointing at ATCA's first re- 
ported coordinates. The observations are compatible with 
SN1998bw lying at that point, and reveal six other objects 
in a (1. X 1. 0) field centered there. As a result of our CB 
model fit to the radio data, as we have explained, we ex- 
pect the optical observations to correspond to SN1998bw, 
and there it is, at the field's center. We also expect, as in 
Fig. ([54I), the CB to be more than two orders of magnitude 
fainter: not observable. It would be nice if this conclusion 
was wrong, that is, if the large "naive" extrapolation from 
radio to optical frequencies in Fig. (^J) was an underes- 
timate by a considerable factor, which is the case for the 
larger extrapolation from optical to X-ray frequencies in 
Fig. (I55I ) (the "naive" prediction there is the one labelled 
p = 2.2). In that case, it may be that a subsequent obser- 
vation of the same field reveals that one of the closer-by 
extra sources has faded away! Three of these sources are 
~ 0. 5 away from the SN, if one of them is the CB, and 
it is dimming, we would not excessively mind that this 
is ~ 60% more distant than the prediction in Fig. (p6|), 
based on a constant-density approximation for the ISM. 

10. The normalization of afterglows 

The values of the CB model parameters that were fitted to 
the broad band data on GRBs with known redshift — and 
are listed in Table III — are narrowly distributed around 
their reference values, except for the overall normalization 
which is much smaller than originally anticipated if 77 = 1. 
This normalization "problem" may point to inaccuracies 
in the various hypothesis that we have made. One exam- 



ple is the contention that all of the incident energy of the 
ISM electrons is radiated in the CB. It may well be that 
the moving CB deflects and scatters the ISM electrons be- 
fore they radiate a large fraction of their acquired energy, 
as suggested by the results of the numerical simulations 
of Achterberg et al. (2001). But abandoning some of our 
simplifications would be premature. Indeed, in conclud- 
ing that in the CB model the normalization of the AGs is 
more than one order of magnitude too large, we have used 
our reference parameters to compute the expected values. 
And there is sufficient elasticity in these parameters to 
obtain a consistent overall picture of all GRB properties, 
as we proceed to review and discuss. 



11. The parameters of the CB model 

With the current analysis of GRB radio AGs we have com- 
pleted a first round of the study of GRB properties in the 
CB model, and it behooves us to look back at the various 
constraints on the relevant parameters. 

In Dar and De Riijula (2000a) we followed Dar (1998) 
and Dar and Plaga (1999) in suggesting that the large 
peculiar velocities of neutron stars may be due to a 
"natal kick" , induced by a momentum imbalance in the 
oppositely-directed jets of CBs accompanying their birth. 
On this basis we chose as a reference value Ecb ^ 10^^ 
erg, for a jet with a reference number of CBs (or promi- 
nent peaks in the GRB light-curve) n^^ ~ 10. Based on 
a first analysis of AG properties, and of GRB 7-ray flu- 
ences and individual 7-ray energies, we set 70 = 10'^ as a 
reference value. 

In Dar and De Riijula (2000b) we investigated two ex- 
treme models meant to bracket the behaviour of a CB as 
it crosses a SN shell, is heated by the collision with its 
constituents, and emits observable 7-rays as it reaches the 
shell's transparent outskirts with a radius R^'j^ , propor- 
tional to its early transverse expansion velocity /3trans c, 
which we assumed to be close to the sound speed in a rel- 
ativistic plasma, /3trans — l/VS- In our "surface" model, 
which is no doubt closer to a realistic description, the en- 
ergy of the GRB in 7-rays is proportional to n^^ [R-cb]^ 7o- 
(Eq. (45) of Dar and De Riijula 2000b). For the chosen 
reference parameters, in the surface model, this predic- 
tion overestimates the GRB fluences by about one order 
of magnitude. Since the individual 7-ray energies corrob- 
orate the choice 70 ^ lO'^, this means that [Rj^'^]^ (and 
Arans) Eire Overestimated by roughly one order of magni- 
tudcfj 

In Dar and De Riijula (2001) we analyzed the X-ray 
"Fe" lines observed in the AGs of some GRBs, which we 
attributed to hydrogen recombination in the CBs, with 
the corresponding Lyman-a lines boosted by a large fac- 



The reduction of thermal free-electron surface-density dis- 
cussed in Sections 6.2 and 6.3 is logically independent from 
the modifications discussed here (such as a decrease in radius) 
which have to do with the fluence produced by the small frac- 
tion of energetic radiating electrons. 
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tor (5/(1 + z) ~ 500. We equated the total number of pho- 
tons in the hnes to the baryon number of the jet of CBs, 
and found agreement with the baryon number in the jet, 
^CB -NcBi to within one order of magnitude. But in the 
current investigation, we have found that the absorption 
of radio waves keeps the CBs hot and ionized (Appendices 
III and V) . This means that our reference value of Nqb is 
likely to be an overestimate. 

In DDD 2001 we proposed a mechanism that would 
quench the expansion of a CB in minutes of observer's 
time, well after it has exited the SN shell. The CBs reach 
an asymptotic radius (Eq. (16) of DDD 2001): 



R„ 



3 NcB /?trans 

2 7rnp 7^ 



(32) 



On the basis of this calculated radius (for /3trans — 
l/(3-\/3), we found that the normalization of optical and 
X-ray AGs agreed with the reference-value expectations. 
On the same basis, we find now that the normalization is 
overestimated by an order of magnitude. The reason for 
the discrepancy is that, in DDD 2001, we effectively placed 
the spectral discontinuity at a "cooling break" frequency 
corresponding to an electron Lorentz factor 70 ~ 1, while 
we have now argued that the discontinuity should occur 
at a higher value 7e = 7(t). 

Both the GRB fluence and the AG fluence are, in the 
CB-model, F oc n^^ , with R = R^^^ for the 7 rays and 
R = Rmax for the AG. At a value of Xqo fixed by the fit to 
the AG's temporal behaviour, the AG fluence is: 



Fag « n^B iic R" 



N, 



1 



CB 



(33) 



All of the above "problems" are solved if we reduce 
our "typical" values of R^ax ^^'^ ^cb by about one order 
of magnitude, relative to our reference parameters with, 
according to Eq. (^), the corresponding reduction of the 
choice of /3trans by half an order of magnitude. 

The precise location of the injection bend is not pre- 
dictable and a modification by up to one order of mag- 
nitude of its position has a small effect on the quality of 
the fits to observations. An increase of the cooling break 
frequency t'b implies a corresponding decrease in AG flux, 
see Eq. (4), adding to the uncertainty in the prediction of 
the precise overall normalization. 

To summarize, the CB model correctly describes, in 
terms of a very limited set of parameters, the properties 
of GRBs and their AGs, including their normalizations. 
This is the case even if we adhere to all of the detailed 
assumptions we have made, even though they are approx- 
imations to a no doubt fairly convoluted physical problem. 

12. Hyperluminal CBs and radio scintillations 

The radio AGs of GRBs often show temporal variations 
of a factor of two or more on a time scale of hours at early 
times and on a time scale of days at later times; e.g. GRB 
000926 in Figs. (|l|, GRB 991216 in Fig. (|l|), GRB 



991208 in Figs. (||, GRB 000418 in Fig. (|5|), GRB 
990123 in Fig. (|j and GRB 970508 in Figs. (|§ l^)- 

Similar variations have not been seen in the optical 
and X-ray bands. The intensity variations of GRB radio 
AGs are very reminiscent of the ones seen in radio signals 
from pulsars in our galaxy, interpreted as scintillations due 
to the motion of the line of sight through the refractive, 
diffractive and dispersive ISM of the Galaxy (see, e.g., 
Lyne and Smith 1982). Some very compact active galac- 
tic nuclei also show an intraday variability that has been 
the subject of much debate (e.g., Wagner and Witzel 1995 
and references therein). At least in one case — the varia- 
tions in the radio intensity of the quasar J1819-I-3845, the 
most extremely variable AGN known at radio wavelengths 
(Dennett-Thorpe and de Bruyn 2000) — it was shown un- 
ambiguously that the variations are scintillations caused 
by the ISM (Dennet-Thorpe and Bruyn 2002). 

The (de)coherence properties in time and frequency 
of the radio scintillations have been used to measure the 
transverse speed of pulsars (e.g., Lyne and Smith 1982). 
Gupta (1995) has demonstrated for a sample of 59 pul- 
sars that their transverse speed, Viss, measured from their 
inter-stellar scintillations, agrees well with their transverse 
speed, Vpni, measured from their proper motion (see also 
Nicastro et al. 2001). 

The movement of the line of sight to pulsars is in most 
cases dominated by their proper motion at a transverse ve- 
locity Vpm larger than the turbulent speeds in the ISM, or 
of the sun relative to the ISM, or of the Earth around the 
sun. The mean Vpm of Gupta's 59 pulsars is 311 kms~^ 
and their mean distance is estimated to be 1.96 kpc. Their 
angular speeds are within an order of magnitude of a cen- 
tral value: 



(Vpm) 

(D) 



5.1 X 10" 



rads 



(34) 



Travelling with a characteristic 7 ~ 10'^ and viewed at 
typical angles of milliradians, CBs have apparent super- 



luminal velocities, V^b of Eq. (|3C|), that are so high (a few 
hundred times the speed of light) that they deserve to be 
called hyperluminal. The angular speed in the sky is: 



>(t) 



VcbW ^, 7(t)<5(t)c( 



Da 



(1 + z) Da 



(35) 



For the reference values 6^1 mrad and an initial 70 ~ 
60 ^ 10"^, the initial angular speed of a CB at redshift 
z = 1 is [I;(^b(0) ^ 2.7 x 10^^^ rad s~^, in the very same 
range as that of Galactic pulsars. The CBs' angular ve- 
locity tj^B (t) and the resulting (inverse) coherence time 
of the scintillations should decline as 7(1) 5{t) does. Both 
pulsars and CBs are pointlikc from the point of view of 
their radio scintillations. Thus, all conditions are met to 
expect pulsar-like scintillations in the radio signals from 
CBs. 

The deviations from a smooth behaviour of the radio 
signals in the case of GRB 980425, as can be seen in Figs. 
( pO| , ^ ), are chromatic, but correlated in time over a much 
longer period than for the other GRBs. Because GRB 
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980425 is so close (z — 0.0085) and is viewed at the unusu- 
ally large angle of - 8 mrad (Table III and DDD 2001), its 
apparent angular velocity, Eq. ( |35| ) , is much larger than for 
other GRBs. The line of sight to this GRB swept a much 
bigger region of galactic ISM than for other GRBs or, for 
that matter, pulsars. Thus, we have no independent infor- 
mation on the ISM irregularities causing scintillations on 
this large scale. 

The analysis of CB scintillations could result in a mea- 
surement of their hyperluminal speeds and a decisive test 
of the cannonball model (fireballs do not have relativis- 
tic proper motions, firecones stop moving close to their 
progenitors and, unlike CBs, have an increasing size that 
should rapidly quench their scintillations with time) . The 
presently available information (or the current informa- 
tion in its published form) is insufficient for us to attempt 
at the moment to extract conclusions about CB hyperlu- 
minal velocities from the observed scintillation patterns in 
the radio AG of GRBs. 

13. Conclusions 

The Cannonball model gives an excellent and extremely 
simple description of all measured properties of GRB af- 
terglows, including their radio afterglows. 

In the CB model, there is an injection bend in the 
spectrum, at the predicted time-dependent frequency z^b 
of Eq. (p^). We have shown that the evidence for the cor- 
rectness of this prediction is very strong, see Table II and 
Figs. (^ and (Q). It is this spectral bend that governs the 
relative normalization of the radio and the optical AG, 
again in agreement with observation, as shown in all our 
figures of AG wide-band spectra. 

Since we have always set the electron index to its the- 
oretical value, p — 2.2, just three intrinsic parameters are 
needed to describe an optical or X-ray AG: 70, Xqo and 
the normalization; the viewing angle 0, although it must 
also be fit, is external to the GRB, like the redshift and 
the absorption in the host and in the Galaxy are. We have 
shown that, in the CB model, the extension of these re- 
sults to the radio domain requires the introduction of just 
one extra parameter: the free-free absorption frequency i^a 
of Eq. (p3|), and that, in spite of various approximations, 
this simplest of descriptions is at the moment entirely sat- 
isfactory. Notice that what one has to parametrize is a two 
dimensional surface: the fluence as a function of frequency 
and time. The shape of this surface is that of a relatively 
simple "mountain" , various cuts of which at fixed t or v 
are shown if Fig. (^). It would be easy, and it may well 
be misleading, to overparametrize this rather featureless 
surface with more than a few parameters. 

It is instructive to compare, or so Occam would have 
thought, the understanding of wide-band AG spectra in 
the CB model with that in the fireball or firetrumpet mod- 
els. In the latter, the number of intrinsic parameters varies: 
seven [e.g. Berger et al. 2001d], eight [e.g. Yost et al. 2001] 
nine [e.g. Yost et al. 2002] and even thirteen [e.g. Galama 
et al. 2000] . This counting does not include the viewing an- 



gle, since the firetrumpets in these works point precisely 
at the observer]^ Moreover, even before the "break" in 
the time-evolution — a period during which it is not in- 
consistent to use the quasi-spherical self-similar approxi- 
mation of Blandford and McKee (1976) for the expand- 
ing material — the ordering of the "breaks" in frequency 
implies a multiple choice of spectral shapes and of their 
evolution (Granot and Sari 2002). 

Countrary to established custom, we are not present- 
ing the values of our fits, which are generally reason- 
able and would become quite good if, again following the 
consuetudinary path, we artificially increased the errors 
to compensate for scintillations in the radio data and/or 
uncertainties in attenuation. The reason is that the CB 
model is a very simplified description of a no doubt very 
complicated reality (e.g. CBs could be somewhat comet- 
like, as opposed to spherical, their inner distributions of 
density, ionization, magnetic field and temperature could 
be non-trivial, even chaotic, etc, etc). Even when the 
physics is much simpler than in the analysis of radio emis- 
sions, and the fits are very good — as is the case in our de- 
scription of optical and X-ray AGs in DDD 2001 — we do 
not report their quite impressive valuesQ. We view our 
"fits" as rough descriptions, rather than true fits. Under 
such circumstances, the overintrepretation of a test has 
every chance of being misleading, much more so in models 
containing many more parameters than the CB model. 

For the same reasons, and because of the systematic 
errors in the data, the values of the parameters we ex- 
tract from our fits should not be taken entirely at face 
value, even though the minimization procedure — which 
attributes to the errors a counterfactual purely statistical 
origin — results in tiny 1 a spreads for the fitted parame- 
ters, and in values that are in most cases satisfactory. 

In the radio domain, as in every other aspect, the 
pair SN1998bw/GRB 980425 is particularly fascinating. 
On the basis of this GRB's observed fluence and distance, 
and given the (totally trivial but all important) depen- 
dence of the fluence on observation angle, we claimed in 
Dar and De Riijula (2000a) that the only peculiarity of 
this pair was that it was observed uncharacteristically far 
from its axis (for a GRB) and uncharacteristically close to 
it (for a SN). In DDD 2001 we proved that the X-ray AG of 
GRB 980245 was also what is expected in the CB model, 
depriving the supernova of its X-ray peculiarity: it did not 
make the observed X-rays. In this paper, by understand- 
ing the magnitude, time- and frequency-dependence of the 
pair's radio signals — which were not emitted, either, by 
the SN — we have demonstrated that SN1998bw was also 

This serious limitation (DDD 2001) is beginning to be 
remedied in the firecone literature, in which the CB-model's 
geometry is being — though with no reference — "standard- 
ized" and its consequences (Dar and De Riijula 2000a) ex- 
plained again; see, e.g. Rossi, et al. (2001), Zhang and Meszaros 
(2001), Salmonson and Galama (2001), Granot J. et al., (2002), 
Panaitescu and Kumar (2001). 

A "x-by-eye" of the figures reporting here the optical AG 
light curves should suffice to prove this statement. 
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"radio normal". Neither this GRB, nor its SN have — in 
the CB model — anything in particular, except the chance 
occurrences of the distance and observation angle. Alas, 
the unique occasion to make a fundamental discovery by 
actually resolving the SN and the CB, as proposed in Dar 
and De Riijula (2000a), may now be very difficult, but, as 
we have explained, not entirely out of the question. 

By pure coincidence, the apparent angular velocities 
of galactic pulsars and cosmological cannonballs are of the 
same order of magnitude. The analysis of radio scintilla- 
tions, one of the methods used to measure pulsars veloc- 
ities, should also be applicable to the GRB ejecta. Thus, 
it ought to be possible to test the CB-model's prediction 
of hyperluminal cannonball velocities. 

Appendix I: The slowdown of a CB 

We review the functional form of the time dependent 
Lorentz factor 7(t), which is explicit and analytical in a 
fair approximation (DDD 2001). 

In minutes of observer's time, CBs reach a roughly con- 
stant radius Rmax and are parsecs away from their progen- 
itor star, a domain where a constant-density ISM may be a 
reasonable approximation. Rclativistic energy-momentum 
conservation in the progenitor's rest frame results in the 
equation governing the deceleration of a CB in the ISM: 

'2 ^72cdtsN- 



9 

, c n 



ttR^ 

P " -^'■max 



(36) 



Interestingly, the above expression is correct both if the 
incoming ISM protons are isotropically reemitted in the 
CB rest frame, or if they are ingurgitated by the CB (in the 
first case, they are reemitted with average energy mp c^ 7^ 
in the progenitor's frame, in the second, the change in 7 
per added proton is d7 = — [mp/M(,3]7). 

Use the relation dtsN = 7(t) dt between the times mea- 
sured in the supernova and CB rest frames, divide both 



sides of the Eq. ( |36D by M^^ 7'^ and integrate to obtain 
the relation: 

1 1 2ct 



7^(t) 



70^ 



where t is CB time, and: 



(37) 



(38) 



with NcB ~ M(,g/mp the CB's baryon number. 

It is important to know the number of electrons accu- 
mulated by a CB as its Lorentz factor decreases from 70 
to 7(t) (in the approximation Ue = Up of a Hydrogenic 
ISM this number equals that of scattered or incorporated 
protons). The number rate of accumulation is related to 
the energy-loss rate of Eq. (p6|) so that: 



dN 



d7 



d7 



(39) 



Assuming constant 77, the total number of ISM electrons 
accumulated at a CB time t is then; 



N(t) =r,N, 



1 



1 

70 



(40) 



The time-dependence of 7(t) with t the observer's time 
is more complicated than Eq. (^^: the relation between 
the two times (dtobs = dtcB(l + z)/(5) introduces a 7 (or 
t) dependence via S. The result for 7(tobs) is (DDD 2001): 



7 = 7(7o,6',Xoo;t) 



C EE 
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C( 



B2 



7o 70 



-bVb2" 

6ct 
(l + z)> 



4616 



1/3 



(41) 



The Lorentz factor of the CB decreases from 70 to 7o/2 as 
the CB travels a distance Xoo/70, whose reference value is 
1.3 kpc, as in Table I. 



Appendix II: The synchrotron cooling break 

We argue that the conventional synchrotron spectral break 
occurs at a very non-relativistic electron energy. The cor- 
responding break in the radio spectrum is unobservable. 

In writing the electron energy-loss rate Eq. (|lO|), 
we have assumed that synchrotron radiation, which is 
quadratic in energy, dominates (inverse Compton scatter- 
ing has the same energy dependence, but it is negligible, 
since the magnetic energy density within a CB is much 
higher than the radiation energy density). The general re- 
sult for the energy loss of high-energy electrons is of the 
form: 



^ ~ Ac (In ^ 
dt mo c^ 



+ a) Ab E + As E^ 



(42) 



The term proportional to Ac describes Coulomb scatter- 
ing and ionization losses in the CB, which are negligible at 
high energies. The second term represents bremsstrahlung 
and its coefficient is (e.g. Shu 1991 ) 



34.35 a ctt c Ub 
2^ 



(43) 



where a = 1/137 and hb is the baryon number density in 
the CB. Adiabatic losses and electron escape would have 
the same energy dependence as bremsstrahlung, but in 
the AG regime we are discussing the CBs are no longer 
expanding and the electron's Larmor radii are so small 
— relative to the CB's radius — that escape losses should 
also be negligible. 

The spectral index of high energy electrons injected 
with a power-law spectrum steepens by one unit at a 
"cooling break" energy [/3(Ec)]^Ec = Ab/As. For 7 ~ 
10'^ and the reference values of fib = lic £ind Uc, the 
synchrotron cooling break is at a subrelativistic energy 
(/3 ~ 0.8). This is in contrast with the injection bend 
at the highly rclativistic energy Eb — 7 mo c^ . The syn- 
chrotron radiation of electrons below the cooling break is, 
for the current data, at unobservably low observer's radio 
frequencies. 
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Appendix III: Photoionization and recombination 
in the CB 

We argue that the synchrotron radiation in a CB is intense 
enough to maintain its plasma partially ionized with ion 
and free electron densities proportional to 7(t). 

The bound-free cross section for photoionization of 
atomic hydrogen in its n-th excited state by photons 
with frequency above the ionization threshold, ;/„ = 
3.29 X 10^^/n^Hz, is given by a,^{n) = ucti gn(iy/t'n)^^, 
with (7i = GAanayiSVS) ~ 7.91 x 10~^^ cm^ (ao = 
0.53 X 10^^ cm is the Bohr radius and gn is the Gaunt 
factor for photoabsorption by hydrogen). For the surface 
flux of photons of Eq. (|^), we obtain an ionization rate of 
the n-th level of atomic hydrogen: 



Ri(n) 



T] Up mp c^7(t)^ n (71 (p - 2) 
8(p + 6) (p- l)hz/n 



(p-2)/2 

■ (44) 
7o, the ionization 



1.1 X 10-^ris-\ 



For the reference values of Uo , p and 7 
rate of the ground state is Ri(l) 

The recombination rate per unit volume of hydrogen 
in an hydrogenic CB is (Osterbrock 1989): 



Rrc 



1.0 X 10"^^ he 



T 



104 K 



-0.7 



(45) 



For the reference value of the CB parameters, and 
T=10,000 K the recombination rate of hydrogen is Rice ^ 
1.1 X 10~^ s~^ . In quasi equilibrium, the ionization rate per 
unit volume, which is proportional to the number density 
of recombincd hydrogen atoms, must be equal to the re- 
combination rate per unit volume, which is proportional 
to the product of the ion and free electron densities in the 
CB. Thus, initially the CB is highly ionized. But, for small 
values of 77 and later times when 7(1) becomes sufficiently 
small, equilibrium between the ionization and recombi- 
nation rates in the partially ionized hydrogenic plasma 
results in lii = lic oc 7(1) . 



For the electron spectrum of Eq. (|13|), the second integral 
in Eq. is dominated by energies Ec ~ Eb (70 7b = 
7). Thus, for a uniform distribution of electrons in the 
CB's volume V, we obtain: 



Hcff 



-^^r(E,) 
V dt ^ ^' 



me 
mn 



377 



Rn 



7- 



(llyycm"'^) 



10' 



n 2 



7 



(48) 



where we have used our reference Rmax- This density could 
be somewhat higher if the emitting electrons are concen- 
trated on the CB's front surface. Note that this density 
increases with time like [7(1)]"^. 

Synchrotron self absorption 

For a power-law distribution, dno/d7o = (p — l)ncff7~P, 
the correct attenuation coefficient for synchrotron self ab- 
sorption at frequency v is (e.g., Shu 1991, Eqs. 19.37,38): 



Xi^ = Kq (p - 1) Ucff B± 



(p+2)/2 j^-(p+4)/2 ^ 



(49) 



where: 
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with ro = e^/mc c^ . An observed frequency I'obs was emit- 
ted at i/(t) = (1 + z)i'ob/<5(t) in the CB's rest frame. 
As an example, z^ob = 5 GHz from a decelerating CB 
with S{t) ~ 7o/2 — 500, if emitted at a typical z = 1, 
corresponds to = 20 MHz in the CB rest frame. For 
our reference parameters, the synchrotron self absorption 
coefficient in the CB is Xi^ — 2-4 x 10^^"rycm~^. Since 
Ucfi oc 7~^, B cx 7, and, after a few observer's days, 
(5 ~ 27 and 7(t) ^ t"^/'^, decreases with time like 



7 



(p+i) ^ t-(P+i)/3 ~ t~i-^ 



Appendix IV: Synchrotron and Bremsstrahlung 
Self Absorption 

We argue that the self-attenuation in the CB of the ob- 
served radio waves is dominated by free-free absorption. 

The density of energetic electrons 

In the rest system of a CB, the ISM electrons arrive at a 
rate d No/dt ~ Up 7 tt R^^x ■ ^ fraction rj of their energy, 
Eb = me c^ 7 , is synchrotron re-radiated. Let be the 
density of the emitting electrons, so that: 

, E. 5^ = / F„ H = / dV dE. E„ ij£ -i^ , («) 

with t(Eo) = Ec/(dEc/dt) the cooling time for electrons 
of energy Ec, determined from Eqs. and ( p^ ) to be: 

r(Ec)=4^ L^. (47) 

mp Up o-^ c 7^ 



Bremsstrahlung self absorption 

The X-ray AG is dominated first by bremsstrahlung from 
plasma electrons and later by synchrotron radiation from 
the swept up high energy electrons (DDD 2001). The ob- 
served X-ray flux, or the theoretical UV flux in the CB 
rest frame, can be used to show that the CB is partially 
ionized during the radio AG observations and that the ion- 
ized fraction of the CB plasma is proportional to 7(1) (see 
Appendix HI). The logarithmic dependence of the plasma 
temperature in the Saha equation on the fractional ioniza- 
tion, keeps the CB's temperature nearly constant during 
its AG phase. Consequently, in the CB rest frame, the 
free-free attenuation at a fixed frequency is proportional 
to [7(t)]^ and the free-free (bremsstrahlung) absorption 
coefficient is that of Eq. ( [20| ) . 

The temperature of the partially ionized CB is of 0(1) 
eV and almost constant during the observed AG. For ^ 
20 MHz emission from a thermal plasma at such tem- 
perature, g '--^ 10 and for one tenth of the typical bulk 
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CB density, lio ^ 10^ cm^'^, one obtains from Eq. ( po| ) 
Xiy — 4: X 10^^" cm^^, which is ~ l-T/yy larger than the 
synchrotron absorption coefficient of the energetic elec- 
trons in the CB (the values of rj are listed in Table III). 
At a fixed observer frequency, = (1 + z) i^oh/S, the free- 
free opacity of the CB decreases roughly like ^ ^"^ ^ t~'^/^ 
compared with the ~ t"^ -'^ decline of the synchrotron self 
opacity. 

The conclusion is that free- free absorption is dominant 
for as long as the ionization of the CB is considerable. 

Appendix V: Emission of self absorbed radiations 

We argue that the energy of the self absorbed radio waves 
and ionizing photons in the CB is radiated mainly by ther- 
mal bremsstrahlung and line emission from the CB con- 
tributing significantly to the observed X-ray afterglow. 

The absorbed radio power is roughly equal to the inte- 
grated emissivity of the CB over all frequencies below z^a, 
defined by Eq. (^). For the spectrum of Eq. and the 
normalization of Eq. (^, the absorbed power is: 

1/2 

• (51) 
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[7(t)]^(p-2) 
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The CB self absorbed radio energy becomes part of the 
thermal energy of the CB plasma. It is radiated by the 
plasma as thermal bremsstrahlung at optical wavelengths 
(in the CB rest frame). For our reference parameters, Va ~ 
100 MHz and rj < 1, this absorbed power is smaller than 
the power absorbed by photoionization. 

The recombination energy is radiated at a rate « 
Rj-ccxNcbl in the CB rest frame where x = ho/iib is the 
fraction of ionized hydrogen in the CB and I — 13.6 eV is 
the binding energy of hydrogen in its ground state. In the 
distant observer frame, the observed radiation is boosted 
and coUimated by the highly relativistic motion of the CB 
and redshifted by the cosmological expansion to: 



dE, 



dt 



R,ecXNcbI(l+z) [6{t)]' 



(52) 



For our reference parameters, hydrogen recombination 
produces X-ray lines with a total energy flux of ^ x^ x 
10~^^ ergs~^ cm~^ . 

Due to their large Doppler shift S, the hydrogen emis- 
sion lines (and the emission lines from the swept up ISM 
and supernova shell material) as well as the CB's thermal 
bremsstrahlung, are shifted to the observer's X-ray band. 
They contribute significantly to the X-ray afterglow and 
may provide a simple alternative explanation (Dar and 
De Riijula 2000) to the commonly assumed Fe-line ori- 
gin of the X-ray lines observed in the afterglows of GRB 
970508: Piro et al. (1998), GRB 970828: Yoshida et al. 
(1999; 2001), GRB 991216: Piro et al. (2000) and GRB 
000214: AntoncUi et al. (2000). 
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Table I - Reference parameters 



Fitted 


Value 


Definition 


e 


10-3 


Observer's viewing angle 


70 


10^ 


Lorentz factor at t = 




1.3 Mpc 


Deceleration parameter 


Other 








10^ 


Dopplcr factor at t = 


Xoo/70 


1.3 kpc 


Distance 70 ^ 7o/2 


R-max 


2.2 X 10^^ cm 


CB's maximum radius 


lie 


10 cm ^ 


CB c number-density 


NcB 


6 X 10^° 


CB's baryon number 


Ambient 






Hp 


10^3 (,ni-a 


Distant p number-density 




lcm-3 


Close-by p number-density 



Comments: The "Fitted" parameters are the typical val- 
ues in the fits to optical and X-ray AGs. The "Other" 
parameters are deduced from the fitted ones (Jo), are cal- 
culated (Rmax), or are deduced from the rest (N^^ and 
he). "Ambient" numbers refer to the ISM, not the CBs. 
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Table III - The Afterglow Parameters 



GRB 


70 


e 


^00 






000301c 


1061 


2.321 


0.128 


552 


0.025 


000926 


787 


0.235 


0.083 


722 


0.027 


991216 


906 


0.403 


0.462 


46 


0.029 


991208 


1034 


0.111 


1.014 


103 


0.011 


000418 


1241 


2.061 


0.332 


298 


0.024 


990123 


1208 


0.464 


0.364 


1604 


0.009 


990510 


1009 


0.261 


0.372 


107 


0.015 


970508 


769 


2.51 


0.516 


559 


0.035 


980425 


495 


7.831 


0.425 


102 


0.007 



Comments: 70: Initial Lorentz factor. 9: Viewing angle 
relative to the CB line of motion, in milliradians. Xo©: 
Deceleration parameter in Mpc (7 = 70/2 at x = Xoo/7o). 
I'a: absorption frequency in MHz in the CB rest frame 
at t=0. T]-. Our best fit normalization divided by the ex- 
pected normalization for the reference parameters that we 
had chosen in previous works. 



Table II - The crossing of the bend frequency 
through the U to I bands [(~ 10 to 3) x 10^^ Hz] 



GRB 


-I 


ti 


/3(ti) 


t2 




970508 


3.7 


0.1-1.5 


-0.58±0.40 


12.1 


-1.12±0.04 


000301c 


5.8 


1.8 


-0.90±0.20 


6-8 


-1.19±0.15 


000926 


7.3 


0.9 


-0.90±0.18 


3.9 


-I.00i0.18 


990712 


13 


0.5-1 


-0.70±0.10 






991208 


17 






3.8 


-1.05±0.05 


010222 


18 


0.20 


-0.88±0.10 


1-5 


-l.lOiO.lO 


991216 


20 


1.67 


-0.58±0.08 






990510 


27 


0.89 


-0.61±0.12 


3.6 


-1.29±0.23 












-1.11±0.12 


990123 


45 


0.033 


-0.69±0.10 


1-3 


-0.90±0.18 



Comments: v^: bend frequency in units of 10^^ Hz. 

/?; spectral index, t;: times after burst in days. For the 
second entry on GRB 990510, see the text. 
References: 

GRB 970508: Galama ct al. 1998a 

GRB 990123: Andersen et al. 1999; Holland et al. 2000 

GRB 990510: Stanek et al. 1999; Holland et al. 2000; 

Beuermann ct al. 1999 

GRB 990712: Sahu, et al. 2000 

GRB 991208: Castro-Tirado et al. 2001 

GRB 991216: Garnavich ct al. 2000a 

GRB 000301c: Jensen et al. 2001; Rhoads & Fruchter 2000 

GRB 000926: Fynbo et al. 2001; Harrison et al 2001 

GRB 010222: Stanek et al. 2001; Masetti et al. 2001 



Table IV - Frequencies, in GHz, at which the radio 

AGs of GRBs of known redshift were measured 



991208 


991216 


000301c 


000418 


000926 




350 


350 










250 






100 


100 


100 




98.48 


86.24 










30 










22.5 




22.5 


22.46 


22.5 


14.97 


15 


15 


15 


15 


8.46 


8.46 


8.46 


8.46 


8.46 


4.86 


4.86 


4.86 


4.86 


4.86 


1.43 


1.43 


1.43 




1.43 


970508 


980425 


990123 


990510 








351 










222 










15 


13.68 




8.46 


8.64 


8.46 


8.65 




4.86 


4.8 


4.88 


4.8 




1.43 


2.49 










1.38 


1.38 







References: 

GRB 970508: Frail et al. 2000a 

GRB 980425: Kulkarni et al. 1998; 

GRB 990123: Kulkarni et al. 1999b; Galama et al. 1999 

GRB 990510: Harrison et al. 1999 

GRB 991208: Galama et al. 2000 

GRB 991216: Frail et al. 2000b 

GRB 000301c: Berger et al 2000 

GRB 000418: Berger et al. 2001a 

GRB 000926: Harrison et al. 2001 
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Fig. 1. — Typical predictions for the CB model's radio 
afterglow. Upper panel: spectra at different times, from 
1 to 300 days. The peak frequencies correspond to CB 
self-opacities of 0(1). The black dots are the location of 
the synchrotron frequency corresponding to the injection 
bend. Lower panel: Light curves at different radio fre- 
quencies, from 350 to 1.43 GHz. The asymptotic curve 
is t~^*^P+^^/'^ (for V 3> as is the case at all frequencies 
shown in this figure). 
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Fig. 2. — Typical predictions for the bend frequency in 
the AG spectrum as a function of time, for 6* = 0, 3/7o and 
10/70. The "optical" U to I band is shown as a horizontal 
band. Upper panel: for 70 — 1250. Lower panel: for 70 = 
750. 
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1 10 

Log time since burst (days) 

Fig. 3. — A comparison between the predicted evolution 
in time of the effective spectral slope through the opti- 
cal/NIR band and the data collected by Galama et al. 
{1998a) for the U, B, V, Rc Ic band of the AG of GRB 
970508 (upper panel), for the K and Rc band (full squares, 
lower panel. Chary et al. 1998) and for the H and Rc band 
(triangle, lower panel, Pian et al. 1998) The three coloured 
lines, in the same order, are the (parameter-less) predic- 
tions. 
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Fig. 5. — Results of a fit to radio and optical observations 
of the light curves of GRB 000301c. The narrowly spaced 
lines refer — from top to bottom — to the K, J, I, R, V, B 
and U bands. The more widely spaced lines refer — from 
top to bottom at the figure's left side — to frequencies 
of 1.43, 4.86, 8.46, 15, 22.5, 100, 250 and 350 GHz. The 
comparison with data is shown in Figs. (^) to (p^). 



Fig. 4. — Comparison between the observations and the 
(parameter-less) prediction for the spectral shape of the 
optical AG of GRB 000301c, at ^ 3 days after burst. Data 
from Jensen et al. (2001). 



S. Dado, A. Dar & A. De Riijula: Radio Afterglows of GRBs 



23 



10 'r 
10' 

10^ 
10^ 
10 

1 

-1 

10 b""' 
10 
10 
10 




1 10 10 

TIME AFTER BURST (DAYS) 



1 10 10 10 

TIME AFTER BURST (DAYS) 



Fig. 6. — Comparisons between our fitted CB model 
AG of GRB 000301c, at z = 2.033, Eq. (|) and Fig. (|), 
with the observed optical data. The figure shows (from top 
to bottom) 1000 times the K-band results, 100 times the 
J-band, 10 times the I-band, the R-band, 1/10 of the V- 
band, 1/100 of the B-band and 1/1000 of the U-band. The 
contributions of the underlying galaxy and an expected 
(but, in this case, unobservable) SN1998bw-Hke SN have 
been subtracted. 
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Fig. 7. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 000301c. Upper panel: the light curve at 350 GHz. 
Lower panel: the light curve at 250 GHz. 
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Fig. 8. — Comparisons between our fitted CB model Fig. 9. — Comparisons between our fitted CB model 

afterglow, Eq. (^), and the observed radio afterglow of afterglow, Eq. (^), and the observed radio afterglow of 

GRB 000301c. Upper panel: the light curve at 100 GHz. GRB 000301c. Upper panel: the hght curve at 15 GHz. 

Lower panel: the light curve at 22.5 GHz. Lower panel: the light curve at 8.46 GHz. 
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Fig. 10. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 000301c. Upper panel: the light curve at 4.86 GHz. 
Lower panel: the light curve at 1.43 GHz. 
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Fig. 11.— The spectrum of the AG of GRB GGG3Glc 
from radio to optical frequencies. Upper panel: in the time 
interval between 1 and 5 days after burst. Lower panel: in 
the time interval between 5 and 10 days after burst. The 
highest peaking curve in the upper pannel corresponds to 
the later time and in the lower panel to the earlier time. 
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Fig. 12. The spectrum of the AG of GRB 000301c 
from radio to optical frequencies. Upper panel: in the time 
interval between 10 and 20 days after burst. Lower panel: 
in the time interval between 20 and 30 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 
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Fig. 13. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed optical afterglow of 
GRB 000926 at z = 2.037. The figure shows (from top 
to bottom) 1000 times the K-band results, 100 times the 
J-band, 10 times the Tband, the R-band, 1/10 of the V- 
band, 1/100 of the B-band and 1/1000 of the U-band. The 
contributions of the underlying galaxy and an expected 
(but, in this case, unobservable) SN1998bw-like SN have 
been subtracted. 
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Fig. 14. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 000926. Upper panel: the light curve at 98.48 GHz. 
Lower panel: the light curve at 22.5 GHz. 



Fig. 15. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 000926. Upper panel: the light curve at 15 GHz. 
Lower panel: the light curve at 8.46 GHz. 
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Fig. 16. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 000926. Upper panel: the light curve at 4.86 GHz. 
Lower panel: the light curve at 1.43 GHz. 



Fig. 17.— The spectrum of the AG of GRB 000926 
from radio to optical frequencies. Upper panel: in the time 
interval between 0.8 and 1.4 days after burst. Lower panel: 
in the time interval between 1.4 and 2.5 days after burst. 
In both cases the highest peaking curve corresponds to 
the earlier time. 
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Fig. 18.— The spectrum of the AG of GRB 000926 
from radio to optical frequencies. Upper panel: in the time 
interval between 2.5 and 4 days after burst. Lower panel: 
in the time interval between 4 and 8 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 



Fig. 19.— The spectrum of the AG of GRB 000926 
from radio to optical frequencies. Upper panel: in the time 
interval between 8 and 20 days after burst. Lower panel: 
in the time interval between 20 and 100 days after burst. 
In both cases the highest peaking curve corresponds to 
the earlier time. 
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Fig. 20. — Comparisons between our fitted CB model 
afterglow, Eq. (||), and the observed optical afterglow of 
GRB 991216 at z = 1.02. The figure shows (from top 
to bottom) 1000 times the K-band results, 100 times the 
J-band, 10 times the I-band, the R-band, 1/10 of the V- 
band, 1/100 of the B-band and 1/1000 of the U-band. The 
contributions of the underlying galaxy and the expected 
SN1998bw-like SN have been subtracted. In a CB-model 
fit, there is in this case some evidence for such a SN (DDD 
2001). 
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Fig. 21. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 991216. Upper panel: the light curve at 15 GHz. 
Lower panel: the light curve at 8.46 GHz. 



Fig. 22. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 991216. Upper panel: the hght curve at 4.86 GHz. 
Lower panel: the light curve at 1.43 GHz. 
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Fig. 23.— The spectrum of the AG of GRB 991216 

from radio to optical frequencies. Upper panel: in the time 
interval between 0.44 and 2 days after burst. Lower panel: 
in the time interval between 2 and 4 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 



Fig. 24.— The spectrum of the AG of GRB 991216 
from radio to optical frequencies. Upper panel: in the time 
interval between 4 and 7 days after burst. Lower panel: in 
the time interval between 7 and 13 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 
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Fig. 25.— The spectrum of the AG of GRB 991216 
from radio to optical frequencies. Upper panel: in the time 
interval between 13 and 30 days after burst. Lower panel: 
in the time interval between 30 and 80 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 



Fig. 26. — Comparisons between our fitted CB model AG 
of GRB 991208, at z = 0.706, Eq. (|), with the observed 
optical data. The figure shows (from top to bottom) 1000 
times the K-band results, 100 times the J-band, 10 times 
the I-band, the R-band, 1/10 of the V-band, 1/100 of the 
B-band and 1/1000 of the U-band. The contribution of 
the underlying galaxy and associated supernova has been 
subtracted. The contributions of the underlying galaxy 
and the expected SN1998bw-like SN have been subtracted. 
In a CB- model fit, there is in this case strong evidence for 
such a SN (DDD 2001). 
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Fig. 27. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 991208. Upper panel: the light curve at 100 GHz. 
Lower panel: the light curve at 86.2 GHz. 



Fig. 28. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 991208. Upper panel: the light curve at 30 GHz. 
Lower panel: the light curve at 22.5 GHz. 
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Fig. 29. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 991208. Upper panel: the light curve at 15 GHz. 
Lower panel: the light curve at 8.46 GHz. 



Fig. 30. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 991208. Upper panel: the hght curve at 4.86 GHz. 
Lower panel: the light curve at 1.43 GHz. 
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Fig. 32. The spectrum of the AG of GRB 991208 from 
radio to optical frequencies in the time interval between 
10 and 14.3 days. The highest peaking curve corresponds 
to the earlier time. 
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Fig. 31.— The spectrum of the AG of GRB 991208 
from radio to optical frequencies. Upper panel: in the time 
interval between 2 and 5 days after burst. Lower panel: in 
the time interval between 5 and 10 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 
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Fig. 33. — Comparisons between our fitted CB model AG 
of GRB 000418, at z = 1.118, Eq. (|), with the observed 
optical data. The figure shows (from top to bottom) 1000 
times the K-band results, 100 times the J-band, 10 times 
the I-band, the R-band, 1/10 of the V-band, 1/100 of the 
B-band and 1/1000 of the U-band. The contributions of 
the underlying galaxy and the expected SN1998bw-like SN 
have been subtracted. In a CB-model fit, there is in this 
case strong evidence for such a SN (DDD 2001). 
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Fig. 34. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 000418. Upper panel: the fight curve at 22.46 GHz. 
Lower panel: the ligfit curve at 15 GHz. 
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Fig. 35. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 000418. Upper panel: the light curve at 8.46 GHz. 
Lower panel: the light curve at 4.86 GHz. 



Fig. 36.— The spectrum of the AG of GRB 000418 
from radio to optical frequencies. Upper panel: in the time 
interval between 9.5 and 30 days after burst. Lower panel: 
in the time interval between 30 and 100 days after burst. 
In both cases the highest peaking curve corresponds to 
the earlier time. 
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Fig. 37. — Comparisons between our fitted CB model 
AG of GRB 990510, at z = 1.619, Eq. (|), with the ob- 
served optical data. The figure shows (from top to bottom) 
1000 times the K-band results, 100 times the J-band, 10 
times the I-band, the R-band, 1/10 of the V-band, 1/100 
of the B-band and 1/1000 of the U-band. The contribution 
of the underlying galaxy and the (in this case unobserv- 
able) associated supernova has been subtracted. 
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Fig. 38. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 990510. Upper panel: the hght curve at 13.68 GHz. 
Lower panel: the light curve at 8.6-8.7 GHz. 
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Fig. 39. — Comparisons between our fitted CB model af- 
terglow, Eq. (H) , and the observed radio afterglow of GRB 
990510. Upper panel: the light curve at 4.86 GHz. Lower 
panel: the spectrum from radio to optical frequencies in 
the time interval between 1 and 6 days after burst. The 
highest peaking curve corresponds to the earlier time. 



Fig. 40.— The spectrum of the AG of GRB 990510 
from radio to optical frequencies. Upper panel: in the time 
interval between 6 and 20 days after burst. Lower panel: 
in the time interval between 20 and 40 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 
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Fig. 41. — Comparisons between our fitted CB model 
AG of GRB 990123, at z = 1.600, Eq. (|), and the ob- 
served optical data. The figure shows (from top to bottom) 
1000 times the K-band results, 100 times the J-band, 10 
times the I-band, the R-band, 1/10 of the V-band, 1/100 
of the B-band and 1/1000 of the U-band. The contribu- 
tions of the underlying galaxy and an expected (but, in 
this case, unobservable) SN1998bw-like SN have been sub- 
tracted. 
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Fig. 42. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 990123. Upper panel: the light curve at 15 GHz. 
Lower panel: the light curve at 8.46 GHz. 
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Fig. 43.— The spectrum of the AG of GRB 990123 
from radio to optical frequencies. Upper panel: in the time 
interval between 0.1 and 1 day after burst. Lower panel: 
in the time interval between 1 and 3 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 
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Fig. 44.— The spectrum of the AG of GRB 990123 
from radio to optical frequencies. Upper panel: in the time 
interval between 3 and 10 days after burst. Lower panel: 
in the time interval between 10 and 20 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 
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Fig. 45. — Comparisons between our fitted CB model AG 
of GRB 970508, at z = 0.835, Eq. (|), with the observed 
optical data. The figure shows (from top to bottom), 10 
times the I-band, the R-band, 1/10 of the V-band and 
1/100 of the B-band. The contributions of the underly- 
ing galaxy and the expected SN1998bw-like SN have been 
subtracted. In a CB-model fit, there is in this case strong 
evidence for such a SN (DDD 2001). 
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Fig. 46. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 970508. Upper panel: the hght curve at 8.46 GHz. 
Lower panel: the light curve at 4.86 GHz. 
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Fig. 47. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 970508. Upper panel: the light curve at 1.43 GHz. 
Lower panel: the spectral behaviour in the time interval 
between 0.12 and 6 days after burst. The highest peaking 
curve corresponds to the earlier time. 



Fig. 48.— The spectrum of the AG of GRB 970508 
from radio to optical frequencies. Upper panel: in the time 
interval between 6 and 20 days after burst. Lower panel: 
in the time interval between 20 and 40 days after burst. In 
both cases the highest peaking curve corresponds to the 
earlier time. 
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Fig. 49. — The spectrum of the AG of GRB 970508 Fig. 50. — Comparisons between our fitted CB model 

from radio to optical frequencies. Upper panel: in the time afterglow of GRB 980425 at z=0.0085, Eq. (^), and its 

interval between 40 and 100 days after burst. Lower panel: observed radio afterglow. Upper panel: the light curve at 

in the time interval between 100 and 470 days after burst. 8.64 GHz. Lower panel: the light curve at 4.80 GHz. 
In both cases the highest peaking curve corresponds to 
the earlier time. 
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Fig. 51. — Comparisons between our fitted CB model 
afterglow, Eq. (^), and the observed radio afterglow of 
GRB 980425. Upper panel: the light curve at 2.49 GHz. 
Lower panel: the light curve at 1.38 GHz. 



Fig. 52. — Comparison between the observed spectrum 
of the AG of GRB 980425 in the radio band and the fitted 
CB model spectrum. Upper panel: in the time interval 
between 2.5 and 20 days after burst. Lower panel: in the 
time interval between 20 and 40 days after burst. The 
highest peaking curve in the upper pannel corresponds to 
the later time and in the lower panel to the earlier time. 
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Fig. 53. — Comparison between the observed spectrum 
of the AG of GRB 980425 in the radio band and the fitted 
CB model spectrum. Upper panel: in the time interval 
between 40 and 70 days after burst. Lower panel: in the 
time interval between 70 and 150 days after burst. In both 
cases the highest peaking curve corresponds to the earlier 
time. 
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Fig. 54.— The late spectrum of GRB 980425 in the 
time interval between 150 and 759 days after burst. The 
highest peaking line corresponds to the earlier time. The 
late isolated point is the last optical observation, which 
must correspond to SN1998bw and not to the AG of the 
CB of GRB 980425. 
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Fig. 55.— The X-ray afterglow of GRB 980425. 
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Fig. 56. — The predicted angular separation of SN1998bw 
and GRB 980425, in milliarcscconds, as a function of time. 



